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PREFACE 

My interest in this subject, interrupted during the war, was 
first stimulated by experimental work which was being carried 
out on solutions of high polymers and, m particular, of t e 
accurate measurements, which were being made at the labora¬ 
tories of the British Rubber Producers’ Research Association, ot 
the thermodynamical properties of the rubber-benzene system 

The theory of solutions of high polymers involves a general 
examination of the statistical mechanics of mixtures, treated as 
co-operative assembUes. While making some contribution to 
the general study of co-operative assemblies, the particular 
application which has been kept in mind has been that to 
solutions of high polymers. It is from this point of view that 
this monograph has been written, although discussion of co¬ 
operative assemblies has been included in so far as it is germane 
to our subject or is necessary in order that the methods which 
have been employed here may be understood. The general aim 
has been to give an account of the mathematical advances in 
the treatment of solutions of high polymers which have been 
made during the last decade. 

The introductory chapter is designed to describe the model of 
a liquid which is used and to recall the nature of a partition 
function and its use in the theory of solutions. The application 
of the particular model to solutions of high polymers and the 
difficulties which arise in the evaluation of their statistical and 
thermodynamical functions is then examined. The next two 
chapters develop the theory and examine its application to the 
available experimental data. Chapter V is devoted to an exami¬ 
nation of the effect of the energy of mixing, and in the final 
chapter various causes to which the remaining discrepancy 
between experiment and theory may be due, and on which 
further research is required, are examined. 

It is a pleasure to acknowledge my indebtedness to Dr. 
E. A. Moelwyn-Hughes and to IVIr. A. B. Pippard, whose kind¬ 
ness in reading the manuscript has ensured greater clarity ol 
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presentation than would othen%ise have been the ca^e. I also 
wish to express my thanks to the officers of the Clarendon 
Press for their unfailing consideration and helpfulness, and my 
appreciation of the skill of their craftsmen. Nor can I cloee thia 
preface \\ithout paying a tribute to the memory of thc^ who 
first guided me in research at Cambridge, and whose deaths 
during the war were not only a personal loss but have left the 
scientific world very much poorer; John Keith Roberts and 
Ralph Howard Fowler. 

A. R. M. 

ROYAli SOCIETi' MO>'D LABORATORY 
TTNIVERSITY OF CAMBRIDGE 

March ms 

ACKNOWLEDGEMENTS 

For permission to reproduce the figures indicated grateful 
acknowledgement is made to the Council of the Cambridge 
Philosophical Society (Figs. 1 and 2); the Council of the Faraday 
Society and the Director of Research of the British Rubber 
Producers' Research Association (Figs. 5, 7, 8, 9, 10 and 12); 
Interscience Publishers Inc. and the Director of Research of the 
British Rubber Producers’ Research Association (Fig. 6). 



CONTENTS 


L INTRODUCTION . ■ • • • • 

Partition functions 
Liquid models 

Partition fimetion of a strictly regular solution 
[I. LONG CHAIN MOLECULES . • • • 

[11. EVALUATION OF THE COMBINATORY FACTOR . 

Bethe’s method ...*•* 
Mixtures containing dimer molecules 
Partial differential equations for the combinatory factor 
Combinatory factor for a mixture containing dimer molecules 
Mixtures containing trimer molecules 

The free energy of mixing , . . • • 

Mixtures containing polymer molecules 

Direct enumeration of distinct configurations 

Basic statistical results . . . . - 

Statistical derivation of general formulae 

Kinetic derivation of combinatory formulae 

[V. COMPARISON OF THEORY WITH EXPERIMENT . 
Vapour-pressure equations ..... 
Measurements on polystyrene .... 
Measurements on rubber ..... 
Discussion of vapour-pressure curves 

The entropy of mixing ..... 

i. THEORY OF MIXTURES WITH NON-ZERO ENERGY 
OF MIXING ...... 

Condition of quasi-chemical equilibrium 

Equivalence of Bethe’s method and the condition of quasi- 
chemical equilibrium ..... 

Free energy of a polymer-solvent system 

Effect of non-zero energy of mixing on agreement with 
experiment ...... 

^I. CRITICISMS AND EXTENSIONS OF THE THEORY . 

The effect of allocating one site each to a benzene molecule 
and to an isoprene unit ..... 

Self-coiling of polymer molecules .... 

The quasi-crystalline model .... 

Statistical methods 


. 1 
3 
7 
11 

18 

28 

29 

33 

37 

39 

39 

44 

45 
48 
52 
55 
59 

G3 

65 

69 

71 

73 

73 

82 

83 

8 () 

88 

97 

1(I3 

103 

106 

no 

111 









I 


IXTRODUCTIOX 

1.1, Liquid mixtures of which one of the components is a high 
polymer behave in a markedly different way from solutions of 
components of low molecular weight. A substance of relati\ ely 
low molecular weight dissolves in a solvent to a precise maximum 
concentration. On the other hand, in the presence of suitable 
liquids, high polymers first swell, imbibing liquid without them¬ 
selves being dispersed. If the polymer is insoluble, nothing 
further happens; in a solvent the high polymer gradually loses 
its form and ultimately disperses to form a solution. For sub¬ 
stances of relatively low molecular weight in the region of dilute 
solutions, various properties such as the vapour pressure, the 
osmotic pressure, and the depression of the freezing-point depend 
in a very simple way on the concentration of the solution which 
is described as ideal. Solutions of high polymers show departures 
from all these simple laws. It is the purpose of this treatise to 
develop the statistical mechanics of solutions of high polymers 
and thus to describe the theory which accounts in large measure 
for these departures from ideal behaviour. First it will be advis¬ 
able to consider briefly the theory of two classes of simpler liquid 
mixtures, usually designated as perfect and regular solutions, of 
which strict definitions will be given below. This procedure has 
the further advantage that it will present in terms of simple and 
familiar systems some of those basic principles of statistical 
mechanics which are to be used in this monograph. 


1.2. The laws of extremely dilute solutions were first eluci¬ 
dated by van ’t Hoff. The equations which he obtained give 
a correct phenomenological description of the relations be¬ 
tween the vapour pressure, the osmotic pressure, and the 
depression of the freezing-point, &c., that is, given any one oi 
these, the deduction of the others is given correctly by van’t 
Hoff's equations. His derivation of these equations cannot be 
regarded, however, as satisfactor\% for his basic assumption is 
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that the solute molecules in dilute solution can be treated as a 
perfect gas, and this is inconsistent -with thermal data.| Tor 
pairs of substances whose behaviour is at least nearly ideal it is 
found that the heat of dissolution of a gas in a hquid is approxi¬ 
mately equal to its heat of condensation; that of a hquid in a 
hquid is zero; that of a sohd in a hquid is equal to its heat of 
fusion. Such results indicate that in solution the solute mole¬ 
cules must be regarded as having a status which is thermo- 
dvnamicaUy indistinguishable from the hquid state. No account 
which starts by assuming that the solute molecules behave as a 
gas can then provide an adequate theory. Furthermore, much 
confusion has arisen as a result of the unjustifiable extrapolation 
of van’t Hofl’s laws to high concentrations. A correct thermo¬ 
dynamical formulation of the laws of ideal solutions was given 
by van Laeri using the methods of Whlard Gibbs. A notable 
step forward was taken when G. N. Lewis§ gave a precise 
definition of a perfect or ideal hquid nnxture and formulated 
the laws obeyed by them. An elementary derivation of these 
laws in terms of thermodynamical cycles was given by Wash- 

bum.il 

It can be shown thermodynamicaUy that necessary conditions 
for two substances to form an ideal mixture are that when 
they are mixed at constant temperature and pressure there is 
no change of volume or energy.tt It is also necessary that the 
different kinds of molecules can be treated as spheres of approxi¬ 
mately the same size; it has been shown that this condition can 
be stated more precisely by saying that the molecular volumes 
of the two species should not differbyafactor of more thantwo.+J 
The behaviour of ideal solutions is summed up most elegantly 

RAniilt/s law. It can be expressed by saying that 


t Compare Kendall, J., Tran^. Faraday 6oc. 66 {lyoi), i. 
t van Laer, J. J., Z. physikal. Chem. 15 (1894), 457. 

§ Lewis, G. X., J- Amer. Chem. Soc. 30 (1908), 668. 

I Washburn, E. W., Z. physikal. Chem. 74 (1910), 385. 
tt It should be noted that zero energj’ change on mixing do^ not req^ 
that the interaction energies be zero. (Cf. Guggenheim, E. A., Proc. oy. 

A 148 (1935), 304, equation (6)). 

Bernal, J. D., quoted by Fowler, R. H., and Guggenheim, E. A.. 

cal Thermodynamics (Cambridge, 1939), p. 351 fm. 
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at constant temperature and external pressure, the partial 
vapour pressure of any component of a liquid mixture bears a 
constant ratio to its molar fraction.f The most satisfactory w ay 
in which to derive these laws is to construct the partition 
function for the assembly and from it to determine the free 
energy. This procedure, which is of quite general application, 
will be illustrated in the succeeding paragraphs. First it is 
desirable to recall the nature of a partition function and to 
describe the liquid model which will be used. 

Partition functions 

1.3. The quantal description of the state of an assembly is 
given completely by the specification of its eigenfunction. In 
the quantal formulation of statistical mechanics, unit weight is 
allotted to each non-degenerate state of the assembly. When we 
speak of a degenerate state of weight m we refer simply to a set 
of states, V7 in number, between which no distinction is made in 
terms of our particular description of the state. If the eigen¬ 
function of the assembly is expressed in terms of the eigen¬ 
functions of its constituent systems, three alternative, and 
familiar, sorts of statistical formulae are obtained according to 
the symmetry restrictions which are imposed on the eigen¬ 
functions of the assembly. On the other hand, if, as is often the 
case, one is concerned only with the macroscopic properties of 
the whole assembly, these alternatives need not be considered 
separately. In describing the equilibrium properties of an 
assembly in thermodynamical terms one is content with a less 
precise description than can be obtained formally from quantum 
mechanics, so that the state of the whole assembly is specified 
by a description of a small number of macroscopic properties. 
In general, at least three such quantities J are required and 
they are to be chosen from the internal energy (E), tempera¬ 
ture (T), pressure (p), volume (F), partial potential (p), and 
number (iV) of particles in the assembly according to which of 

t Henry’s law is a special case of this. 

t Or two quantities if the masses of the components, and therefore the 
parameters N, are fixed. 
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these quantities it is convenient to choose as independent 
variables.! 


1.4. In the most common specification of an assembly the 
temperature, volume, and the numbers of molecules of the com¬ 
ponents are taken as the independent variables. If tcr^ is the 
degeneracy of the quantum state of energy then the sum of 
terms exp(—€,./feT), each we^hted according to the degeneracy 
of the corresponding eigenstate, and summed over all the states 
of the assembly, is the (ordinary) partition function for the 
systems of which the assembly is composed. Thus 

/= 2®-rexp(—(1.1) 
r 


The partition function, which is a dimensionless quantity, sum¬ 
marizes in a convenient mathematical form the way in which 
the energy of the assembly is partitioned among the individual 
systems of which it is composed. It is the analogue of Gibbs’s 
canonical ensemble. The related thermodynamical potential 
is the Helmholtz free energy which is given by 



The corresponding expression in terms of the macrodesciiption 
of the assembly gives a well-known formula 


/27r7HfeT\|p. 

\~hr-j 


for the partition function per molecule. We must now consider 
the partition function for the whole assembly. 


1.5. The energy of an assembly is the sum of (o) the kinetic 
energy of translation of the individual molecules or atoms 
(systems) of which it is composed, (6) the rotational, vibrational, 
and electronic energies of the individual systems, and (c) the 
configurational energj'^ of the assembly as a whole, due to the 
interactions between its component systems. Corresponding 

t There may sometimes be additional parameters required; e.g., the fiactioni 
of a given quantity of hydrogen which are in the para and ortho states, or thi 
fractions of a given quantity of a component which are liquid and vapour. 
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to thisf the partition function F(T) of the assembly splits into 
three factors. Thus 

F{T) = (1-3) 

where the partition function corresponding to the translational 
motion of a single particle is given by the product oU{T) and the 
configurational (volume) factor per system, which is here in¬ 
cluded in j(T) is the partition function of a single particle 
corresponding to the energy which arises from all its internal 
degrees of freedom, andQ(T) is the configurational term for the 
whole assembty. If the rotational and vibrational motions may 
be treated as independent, then j(T) splits into factors: r{T) for 
the rotational energy, v(T) for the \dbrational energy, and e{T) 
for the electronic energy. 

1.6. In the applications which are to be discussed in the 
following chapters, we are to consider the equilibrium of matter 
in bulk. Consequently, we shall be concerned to distinguish 
only between different macrostates of the assembly, lumping 
together all microstates for which the whole assembly has the 
same configurational energ}’. In effect, this assumes that l(T) 
and j(T) are independent of the configuration of the assembly 
as a whole. In other words, there are g( IT) states of the assembly 
with configurational energy W between w'hich we do not dis¬ 
tinguish in terms of the particular description of the assembly 
which is relevant to our purpose. Thus the determination of 
the partition function reduces to the evaluation of the con¬ 
figurational factor Q(T). In somewhat less general terms it can 
be expressed thus: consider an assembly, specified by the para¬ 
meters N which are relatedj to the numbers of molecules of the 
different species. Let the assembly be in one of 5 r(N') configura¬ 
tions in each of which the potential energj^ is ir(A"). The con¬ 
figurational term is then given by the multiple integral 

n(T) = j...j exp{- W{N)}JkT n {d<or, (1.4) 
where n is an abbreviation for the product of elements of 

t See § 1.10, equation (1.9). 

i In addition to the numbers of particles of different kinds, there will also be 
parameters which give the numbers of pairs of neighbours of different kinds. 
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volume for each type of molecule, equal in number to the num¬ 
ber of molecules of that t^^pe, and is used to indicate that the 
integration is to be extended over all the space and velocity co¬ 
ordinates of all the systems. From what has been said already 
(§ 1.4), it is clear that, in terms of a description of the assembly 
relevant to our purpose, the partition function must consist 
formally of the sum of terms of the type 

y(.Y)exp{- ir {X)lkT} (1.5) 

apart from factors arising from the effective free volume and the 
partition functions for the translational and internal degrees of 
freedom of each molecule. The essence of the problem is to 
determine the appropriate combinatory factor g{N). To write 
an expression like (l.o) merely shifts the point of difficulty from 
the evaluation of the multiple integral in (1.4) to that of the 
appropriate combinatory factor g(N). 

1,7. A generalization of the idea of a partition function was 
given by Fowler, f who considered the statistical analogue of 

a very powerful 

method and is usually called the grand partition function. In 
using the grand partition function the temperature, volume, and 
partial potential are taken as the independent variables, and the 
appropriate weight factor in the summation over all the states 
of the assembly is exp{—(e,—A/x)/fe7’}. The grand partition 
function is defined by 

exp{-e,lkT)X^, ( 1 . 6 ) 

N r 

in which the sum is taken over all states of the assembly and over 
aU values of the numbers of systems of each kind composing it. 
In this expression A is the absolute activity which is defined by 


Gibbs’s grand canonical ensemble. It provides 



fi = kTlogX, 



where G is the Gibbs function. The ‘thermodynamical potential’ 


related to the grand partition function is given by 



felogZ = klog 2 2 

N .r 



€,/kT)X^. 


t Fowler, R. H., Proc, Cambridge PhU. Soc, 34 (1938), 382. 
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The corresponding expression for the grand partition function 
in terms of a macrodescription of the assembly is 

g{W)ex^i-WlkTm)^, ( 1 - 8 ) 

N W 

where ^ includes the volume factor and the partition function 
per system corresponding to the translational and internal 

degrees of freedom. 

This brief descriptionf of the nature of a partition function 
has been included because we shall have constant recourse to 
those methods of statistical mechanics, due to Fowler and Dar¬ 
win, which use partition functions. 

1.8, The simplest assemblies of statistical mechanics are 
those in which the interactions between the elementary particles 
of which the assembly is composed are so weak that they can be 
neglected. In this case the assembly can be regarded as a large 
number of virtually independent systems, and it is called an 
ideal or perfect assembly. If the interactions, while no longer 
negligible, are still small, one has a slightly non-ideal assembly. 
However, there are other assemblies in which the interactions 
are so important that the states of any distinct system are 
influenced in a fundamental way by which states of the other 
systems are occupied. Such assemblies are termed co-operative. 
The constituent systems of a co-operative assembly carmot be 
treated as independent, that is, the free energy (for example) 
is not simply the sum of the free energies of the constituent 
systems. The advances which have been made during the last 
decade in the application of statistical mechanics to liquids have 
followed from the recognition that a liquid forms a co-operative 
assembly. 

Liquid models 

1.9. Earlier discussions of liquids involved an unwarranted 
extrapolation of results which had been derived only for slightly 
imperfect gases. Such extrapolation was based on the assump¬ 
tion that a liquid does not differ in any essential way from a gas, 

t For a more detailed description see Fowler, R. H., Statistical Mechanics 
(Cambridge, 1936), ch. ii; Fowler, R. H., Proc. Cambridge Phil. Soc. 34 (1938), 
382; Guggenheim, E. A., J. Chem. Phya. 7 (1939), 103. 
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being merely more highly compressed. It is correct that there 
is no discontinuity between gas and liquid, but such a formal 
identity overlooks the fact that an essential part of the treat- 
mentf of slightly imperfect gases is that the effective volume 
is not very much smaller than the total volume. Even at the 
critical point the effective volume is only about half the actual 
volume and in hquids it becomes even smaller. That is, the 
usual treatment of imperfect gases requires the assumption that 
the difference between the total and effective volumes is very 
much smaller than the total volume, while in a liquid the reverse 
is true, namely, the effective volume is very much smaller than 
the total volume. 


The simplest useful model of a liquid of symmetrical mole¬ 
cules is a collection of spheres packed together, not tightly (which 
would represent a crA’stalline solid), but not very far from tightly. 
It is clear that such a model has less order than a crystalline 
solid, but it retains a degree of order far different from the 
almost complete randomness of configiuation which is character¬ 
istic of a gas. This model will be referred to as quasi-crystaUine. 
Such a model ensures that the component molecules are allowed 
considerable freedom of movement; at the same time it secures 
a high degree of local order. That there is a high degree of local 
order in certain Liquids at ordinary temperatures has been con¬ 
firmed by X-ray analysis and it appears that, except in the 
neighbourhood of the critical point, a liquid is much more like a 
ciystal than a gas.§ At temperatures weU below the critical 
temperature the number of neighbours surrounding each mole¬ 
cule has a well-defined average value. The fiuctuations about 


t With the exception of the treatment of condensed systems due to Mayer 
(Mayer, J. E., J. Chem, Phys. 5 (1937), 67; Mayer, J. E., and Ackerman, P.G., 
ibid. p. 74; Mayer, J. E., and Harrison, T. F., ibid. 6 (1938), 87, 101). 

* Friedrich, W., Phys, Zeit, 14 (1913), 397 ; Debye, P., and Scheirer, P., GoU, 
}iach. 16 (1916); Debye, P., and Menke, H., Phys, Zeit, 31 (1930), 797; Ergeb. 
d. Tech. Pordgtnk. 11 (1931): Bragg, W. L., The Crysialline State, A General 
Survey (G. Bell and Sons, 1933), p. 195; Gingrich, N. S.. Eev, Mod, Phys, 15 
(1943). 90; Finbak, C., Skr. norske Vidensk, Akad, 10 (1945), 

§ This viewpoint has recently been elaborated and emphasized by Frenkel, 

J. . Kinetic Theory of Liquids (Oxford, 1946). See also Meyer, K. H., Zeits.f, 
physik. Chernie, B 44 (1939), 383; Helv, Chim, Acta, 23 (1940), 1063; Meyer, 

K. H.. and Boissonas, C. G., ibid. 23 (1940), 430. 
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this average value are not serious; except for a rather unimpor¬ 
tant fraction of centres in the liquid, they will not disturb the 
regularity of the geometrical relationship between each molecule 
and its closest neighbours, which, on the average, is very similar 
to that in a crystal. The fluctuations are, however, usually 
sufficient to destroy the long-range order, so that there is no 
regular geometrical relationship between a molecule and its 
more distant neighbours. On this model, which is the simplest 
that is in agreement with the available X-ray data on liquids, 
the essential difference between a liquid and a crj’stal is that 
whereas the former possesses a high degree of local order, w hich 
can be defined by an average coordination number, the latter, 
in addition, exhibits long-range order in the geometrical relation¬ 
ships between its constituent molecules. 

1.10. It is a standard and well-known result of statistical 
mechanicst that the partition function of an assembly of N 
molecules of mass m which are confined to a volume V is 

_ ^nrp)j{T)}m(T), (1.9) 

where the configurational term Ci{T) is given by 

Q(T) = ^ |... J (1.10) 

Here j(T) is the partition function for all the internal degrees of 
freedom (rotational, vibrational, electronic, and nuclear spin); 
1{T) is the partition function for the translational motion of a 
molecule; TF(.ft'^) is the configurational potential energy of the 
assembly referred to that of the molecules at infinite separation 
as zero; the factorial is introduced to remove multiple counting of 
indistinguishable configurations; and (dw)^ denotes the product 
of the N elements of volume to which the centres of the N mole¬ 
cules belong. The multiple integration is extended over the whole 
of the 3N dimensional phase space accessible to the molecules. 
For an assembly of non-interacting spherical molecules or for 

t Fowler, R. H., Suaisticol Mechanics (Cambridge, 1936), § 2.6. 
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a liquid, neglecting molecular interactions, it can be shown 
thaty nirr\ __ /,- \v 


a{T) = (evy\ 


( 1 . 11 ) 


where v is equal to the effective free volume per molecule and is 
a function of VIN. This, however, can be of little practical 
interest on accoimt of the neglect of the intermolecular forces. 


1 . 11 . For a liquid, the volume F in which it is contained is 
deter mine d by the intermolecular forces. A useful approxima¬ 
tion to take account of the interactions in a normal liquid is 
provided by the smoothed potential model, in which the con¬ 
figurational energy is replaced by its value averaged over all 
accessible configxu^tions. This average value of the configura¬ 
tional energy can be written as 



where —x is the average value of the molecular energy for all 
relative configurations of the given molecule and its neighbours; 
since the attractions are assumed to be short-range forces, only 
near neighbours will contribute appreciably to it, so that its 
value will be determined by the average number of near neigh¬ 
bours within a specified distance. It therefore follows that x is 
a function onlj’ of the average volume per molecule. This model 
leads to the familiar result 


n{T) = (ev)-''exp(Ax,fer) - [uexp{(x-ffeT)/feI^]-\ (1.12) 


In general, the configurational term Q.{T) given by the integral 
in equation (1.10) will reduce to the product of three factors. 
Thus it will be the product of a combinatory factor which arises 
from the number of distinguishable configurations of the systems 
of the assembly, a volume factor, such as (ev)^\ and a term which 
depends on the configurational energ}', such as exp(A’x/fe?’). 

Derivations of the configiirational term which are based on 
more precise physical models can be formulated. A model, due 
essentially to Mie,J in which the change in interaction energy 
as the molecule moves about its equilibrium position is allowed 
for b}’ assuming that each molecule moves in a field correspond- 


t Fowler, R. H., SlatiMical Mrchanics, 2nd edition (Cambridge, 1936), § 13.1. 
i Mie, G., Ann. d. Physik, 11 (1903), 657. 
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ing to an isotropic three-dimensional harmonic oscillator, pro- 
rides a useful comparison with a solid. A more refined modelj* 
in which the molecular energj^ is calculated directly from the 
interaction energy between pairs of molecules has also been 


examined. The factor e in equations (1.11) and (1.12) represents 
an essential point of difference between a liquid and a solid. It 
gives rise to an additional term k in the molecular entropy and 
specifies the greater disorder of the liquid. J 


Partition function of a strictly regular solution 

1.12. A mixture of molecules of two species i and j is said to 
form a strictly regular solution§ if (i) both types of molecule 
pack in the same way so that they have a common coordination 
number; (ii) the molecular volumes are sufficiently alike for a 
mixture of the two kinds of molecules to pack in the same way as 
each of the pure liquids; (iii) the ratio of the free volumes of the 
two pure liquids lies between 1 and 1*30; (iv) both the molecular 
volumes and the free volumes are unaltered when the tw'o liquids 
are mixed at a given temperature and pressure, and (v) only 
closest neighbour pairs make any significant contribution to the 
potential (configurational) energy of the mixture. 

An energy of mixing is defined so that if, in the two pure 
liquids, an interior molecule of type i is interchanged with an 
interior molecule of tj^e j the total increase in the potential 
energy is 2w^^. If the average potential energy of a molecule of 
tj^e i in the pure liquid is — Xt that of a molecule of type j 

in its pure liquid is — X;> then the average energies of interaction 
of a pair of closest neighbour molecules i is — 2 x ,/2 and of a pair 
of closest neighbour molecules j is —2xjlz, where z is the common 
coordination number of each pure species and of the mixture. 
If the average interaction energy of a molecule i and a molecule 


t Lennard'Jones, J. E., and Devonshire, A. F., Proc. Roy, Soc, A 163 (1937) 
63; ibid. 165 (1938). 1. 

t The discontinuity between the thermodjuiamical functions of liquids and 
crystals was first emphasized by Eyring (Hirschfelder, J., Stevenson, D.. and 
Eyring, H,, J. Chem. Phys. 5 (1937), 897) and by Lennard-Jones (Lennard- 
Jones, J. E., and Devonshire, A. F,, Proc. Roy. Soc. A 163 (1937), 61). 

§ Fowler, R. H., and Guggenheim, E. A., StatMcal Thermodynamics 
(Cambridge, 1939), ch. viii. 
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j which occupy a pair of closest 
are in contact) is Ty^-. then 


neighbour sites (that is, which 
-Tf—(1.13) 


where T),-. Vjj have been wTitten for —'2xi Z and —2xj Z respec- 
tivelv. 


1.13. Consider a particular configuration of the mixture of 
S^ molecules of type i and molecules of t^-pe j in which there 
are zX closest neighbour i — j contacts. The numbers of i—i and 
j—j closest neighbour contacts are then Z)and i) 

respectively. The potential energy is therefore given by 


ir = -X 


(hl4) 

The configurational factor in the partition fimction of the mix¬ 
ture is, from equation (1.10), given by 


0.(T) = r f 

J J 

(1.15) 

First consider the special case of a perfect solution defined by 
zero; then 


Q(T) = 




1 


: . .A.» c ♦ 


(1.16) 


1* 


With an error of not more than about 1 per cent, the arithmetic 
mean of the free volumes may be replaced by the geometric 

mean. Thus v ^ v 

and with this approximation equation (1.16) reduces to 


n(T) 




y , {ft^, exp(x,- kT)y^€Vjex^{xj 

(1 17) 


gi\'ing the configurational factor in the partition function for a 
perfect solution. It should be noted that the first factor in 
equation (LIT) i.s the number of ways of an-anging iN,- indis- 
tinc^uishable svstems of one kind and X indistinguishable 

w ' 
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systems of a second kind on an array of N^+Nj sites. For the 
grand partition function of a perfect solution this gives 

Z = 2 g(N„N,)<>xtl(N,x,+N,x,)lkT}{a,MT)UT)y‘X 

x{£u,yJ’)i,(J’)W'Af<, ( 1 . 18 ) 

where g{Ni,Nj) is the combinatory factor, that is, the number of 
distinct configurations of the systems of which the assembly is 
composed with an energy W equal to 
example of a perfect solution giN^, Nj) is equal simply to 

{N,+N^)W.N^l 

The partition function for a pure liquid is, from equations 
(1.9) and (1.12), 

{l{T)j(T)ex^(xlkT)ev}^ = {L(T)}^. (1.19) 

Then the partition function for a given configuration of 
molecules i and molecules j which form a perfect mixture can 

always be expressed by {Li{T)}^\L^(T)}^'i provided L^iT) and 
Lj{T) are independent of the ratio NJNj. The partition function 
for the system of molecules i and Nj molecules j including 
all configurations is therefore 



{N,+Nj)\ 

Ni\Nj\ 



( 1 . 20 ) 


in agreement with the form obtained (equation (1.17)) by evalua¬ 
ting the integral in equation (1.16). The grand partition function 
for a perfect solution is then given by 



2 


{Ni+Nj)\ 


{Lj{T)Xjy’^<{Lj(T)Xj}^’s. 


( 1 . 21 ] 


1-14. The laws of perfect solutions can be derived immedi¬ 
ately firom the partition function. From equations (1.3) and 
(1.17) the free energy of the mixture is given by 


A^-{log (T) -|-log j,( T)+log u,-} —Ni—Nf XilkT— 
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From equations (1.9) and (1.12) the free energy of the pure 
liquid i is given by 


‘ii 

kT 




so that the free energy of mixing is given bv 


A.4' 

kT 


a;.)+.v, log log .V.. (i.23) 


Denote by /xj the partial potential of the species i in the liquid 
mixture, and by pj „ the partial potential of the pure species i. 
If the Gibbs free energy of mixing is denoted by IG, then for a 
binary' mixture it is defined bv 

G = 

where and G^ are the Gibbs free energies of the pure species 
and G is that of the mixture. At constant pressure and tempera- 

m =a) +(c“) 


Thus 


if T.p 


Mi 


-I 

Mi.OT 


‘^4) (.u,) 


C 


Thus, provided the pressure is small. 


Mi = 

- iJio^kT log ’ 

1 

(1.24) 

Mi = 

-Mlo-r*riogy 

(1.25) 


The free energy of the gas phase is given by 


A 


kT 




yU'vmiK 


Y'. 


4-1 


(1.26) 


where primes are used to specify parameters of the gas phase. 
If the vapour of each species can be treated as a perfect gas, then 
the partial pressures in the gas phase are given by 


Pi 


X\kTi V' (all i). 


(1.27) 
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From equations (1.26) and (1.27) the partial potential in the 
vapour phase is given by 


/4 


kT\og{UT)3\{T)kTlp,), 


from which it follows that 


= Mi,0+*^ log 


( 1 . 28 ) 


Pifi 


where the subscript 0 is used t< 
equilibrium it is necessary that 
phases should be equal. It foll( 

Pi = Pifi 


K 


Ni+N/ 

K 


(1.29) 


Pi 


which is Raoult’s law. Equations (1.29) have been obtained by 
assuming that the energy of mining is zero, without recourse to 
the more restrictive condition of negligible interactions between 
all pairs of molecules. 

1.15. It has been shown that, in general, the configurational 
factor in the partition function for a strictly regular solution is 
given by 

il{T) = J... J 

(1.15) 

The evaluation of the integral in equation (1.15) can be carried 
out only approximately. A quantity X defined by the equation 

e-fu>^kT j_j = J...J e-^“«"^(da,,.)^*(d<.,^)-v, 

(1.30) 

can be introduced. The integral in the left-hand member can be 
evaluated as in equation (1.16); this leads to the result 


il(5r) = jeu^exp 


w 


cu^exp 


tell 




exp 


Xw 

(1.31) 




1(3 


I^'TRODUCTIOX 


-_ Cii. I, 5 1.15 

this was the device used bv Gueeenheim.t althou 2 h he neglected 

^ w w w 

to distinguish between A' defined by equation (1.30) and the 
average or equilibrium value X. However. X and X are not 

identical! and the relation between them can be seen in the 

» 

following way. By clitferentiating equation (1.30) with resi>ect 
to T it can be .<hown that, provided is sensibly independent of 
T. then 'r\ 


X = X-T 


cX 

cT 


c{X T) 

f(l 7)’ 


(1.32) 


where A' is the equilibrium value of A and is defined, in the usual 
wav. bv 

(1.33) 

Alternatively, in the light of what has been said m §§ 1.6 and 1.11, 
QfT) can be expressed formally as 

0(7) = g(.\..A;. A)exp(-Ai<-,j *7){fa;exp(x,- fe7)f‘<x 

X (euj exp(xq kT)}-''!, { 1.34) 


where g( A'. A’. A) is the number of distinguishable configurations 
of the Aj- molecules i and A^ molecules j in wiiich there are rA 

closest neighbour i—j contacts. 

To proceed farther it is necessary, on the one hand, to evaluate 
A or A or. on the other hand, to evaluate A,-^-). Xo 

simple method is available, and it is necessary to fall back 
on approximate methods, the detailed consideration of 5^hich 
is deferred to a later stage (Chaptez's III and ^ ). In fact, 
the difficulty in the evaluation of the configurational term in 
the partition function for a strictly regular solution has led to the 
kernel of the difficulty of the statistical mechanics of solutions 
of hiiih polvmers. For when the assembly contains particles of 
markedly different size, the combinatory factor, even when the 
energy of mixing is zero, contains parameters which involve 
pairs of sites and their interrelation, in a way similar, for example, 
to the parameters A of equations (1.15) and (1.34). Thus the 

T Gus^enlieiin. E. A.. Proc. Roy. Soc. A 148 (1935), 304. 

t Koshbrooke, G. S., ibid. A 166 (193S), -96. 
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determination of the partition function for mixtures containing 
polymer molecules, whether the energy of mixing be zero or not, 
involves problems of the same kind as that imphcit m the 
evaluation of equations (1.15), (1.31), or (1.34). The nature and 
occasion of the essential difficulty will become clearer by con¬ 
sidering solutions of long chain molecules in relation to the quasi¬ 
crystalline model of a liquid which has been adopted. 


I 



II 


LONG CHAIN MOLECIHES 

2.1. The quasi-crystaUiiie model of a liquid which has been 
described in ^ 1.9 will be used. On the averagfe there is a definite 
relationship between any molecule and its near neighbours 
which can be specified by a coordination number although the 
fluctuations about the average positions are sufficient to destrov 
any long-range order. From another point of view the whole 
hquid can be regarded as providing a potential field of which the 
minima provide an array of sites which can be occupied by the 
molecules of different kinds. There is a definite geometrical 
relationship between any site and its near neighbours which, on 
the average, is ver\- similar to that in a crystaL but there is no 
long-range order. For a mixture in which the molecules of each 
kind are at least roughly spherical and of about the same size, 
each molecule of whatever kind can be regarded as occupying 
one site. The combinatory factor in the partition function arises 
from the number of aixangements of molecules of different kinds 
over the available sites for each value of the interaction ener^. 
This is precisely what has been done in considering regular 
solutions in 1.13 and 1.15. The task now is to examine the 
use of the quasi-crystalline model of a liquid when the molecules 
of one species are long chains. 

2.2. First, consider a simple chain molecule which is formed 
by the repetition of a segment or sub-molecule. Any number of 
such substances can be called to mind; one which occurs naturally 
is rubber, in which the sub-molecule is an isoprene unit. If the 
whole molecule consists of r units (r ^ 1). then a suitable repre¬ 
sentation of it is to consider it to occupy r sites of the array in 
such a way that adjacent sub-molecules in the chain occupy 
closest neighbour sites. Each sub-molecule is regarded as 
occupying one site on the array, and this allows not only for the 
fact that a chain molecule is being considered but also permits 
a hio’h degree of flexibilitv of each chain molecule. It is clear that, 
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CHAIN MOLECULES 


by the aUocation of an appropriate number of sites, in which the 
sites which accommodate any pair of adjacent sub-molecules 
are suitably disposed, both branched chain and closed rmg 
molecules can be discussed in terms of the quasi-crystaUine 
model of a liquid. A physical representation of this sort also 
makes it possible to consider mixtures of species of which the 
volumes of the molecules differ in the ratio of small integers. 
For convenience and brevity, chain molecules* which contain 
2, 3, 4,..., r,... segments or sub-molecules will be referred to as 
dimer, trimer, tetramer..., r-mer,... molecules. In conformity 
with this convention, a simple molecule which can be accommo¬ 
dated on one site of the array will be called a monomer. 

2.3. Consider a molecule which contains r units. The number 
of closest neighbour pairs of sites formed by the sites occupied 
by the r-mer and the rest of the sitesf is denoted by zq. This 
quantity is the number of closest neighbour contacts or inter¬ 
actions into which the polymer molecule enters, for the defim- 
tion excludes pairs of closest neighbour sites which are occupied 
by adjacent units of the polymer. We will now examine the 
value of zq for different kinds of polymer molecules. J For a 
simple chain molecule which contains r groups in aU, there are 
z—2 closest neighbom* interactions for each internal group and 
z—1 for each of the end groups. Thus there are 

(r-2)(z-2)-f2(z-l) = rz-2r-f 2 

closest neighbour interactions in all. For a branched chain 
molecule there is one less interaction at the jimction of each 
branch, but each branch introduces an additional end group 
and so an additional interaction. Thus the total number of 
interactions for a branched chain molecule is the same as for a 


simple chain 
q’ given by 


olecule. In each case the parameter q has a value 


{r-q’)z = 2(r-l). 


( 2 . 1 ) 


For a closed ring molecule this no longer holds. The essential 


t If the molecule coils back on itself so that two non-adjacent units of the 
molecule occupy a pair of closest neighbour sites, this pair should be included 
with regard to each of the sites. 

t Miller, A. R., Froc, Cambridge Phil. Soc. 43 (1947), 422. 
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fact can be expressed by saying that a closed ring molecule has 
no end groups; each group in a closed ring molecule gives rise to 
z—'2 closest neighbour mteractions. Thus for a closed ring mole¬ 
cule the parameter q has a value q" given by 

(r-g'): 2r. ( 2 . 2 ) 

It should be noted that for large r (r > 1 ) the values 9 ' and 9 ' do 
not differ greatly. In the following paragraphs wherever the 
parameter 9 is used it is implied that it is to bt^ given its appro¬ 
priate value 9 ' or 9 ' according to the t}^ of molecule. 

2,4. It will be useful at this stage to derive certain statistical 
relations which will be of frequent use. Consider a mixture which 
contains molecules of species i (all i); each molecule of species 
i contains segments and occupies r,- sites of the array in such a 
way that adjacent segments occupy closest neighbour sites. A 
parameter 9 ,- is introduced for each species and is defined by the 
fact that each molecule i shares in 29 ,- closest neighbour inter¬ 
actions. The precise values of 9 , are determined by equations 
( 2 . 1 ) and ( 2 . 2 ), but are unimportant for the immediate purpose 
and play no part in the deduction of the results of this section. 

Quantities (all i and j for which i ^ j) are defined so that 
zXij is the number of closest neighbour pairs of sites of which 
one is occupied by a molecule of species 1 and the other by a 
molecule of species j. Quantities J,-.- (aU i) can also be introduced 
and defined by the fact that zX,-,- is the number of closest neigh¬ 
bour pairs of sites each of which is occupied by a molecule of 
species i including those closest neighboxir pairs which (due to 
the molecule coiling back on itself) are occupied by non-adjacent 
units of the same molecule of species i but excluding those cases 
in which the pair of sites is occupied by adjacent units of the 
same molecule. A further quantity, Xa, is defined as follows: 
let zx^i (aU i) be the number of pairs of sites which are occupied 
by two adjacent units of the same molecule of species i. 

A molecule of species i shares in two kinds of closest neigh¬ 
bour interactions. It may give an interaction with another 
molecule i or with a non-adjacent sub-molecule of itself (forming 
an 2 —i interaction in each case), or it may give an interaction 

Th9 J 3 K University Library 
Acc. No. ^ ' 
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with a molecule of another species (of which an i-j mteraction 
is typical). In totalling all the interactions in which all the mole¬ 
cules of species i take part each of the former kind will be in- 
J latter kind will be included once. 


Hence 


2^,i+ 2^ ^ik 

k 


q^Nf (all t). 


(2.3) 


Consider the r,- sites which accommodate a molecule of species i. 
Each site has 2 : closest neighbours. The difference between 
and ZQi is simply twice the number of closest neighbour pairs 
which are occupied by adjacent sub-molecules of the same mole- 

cule i. Hence 


2 ' ^ik 
k 




(aU t). 


(2.4) 


2a:,.,. 


From equations (2.3) and (2.4) it follows that 

= (U-qM (all i) (2.5) 

for all types of molecules. It is an easy matter to express the 
quantities a:,, in terms of the alone, but these expressions have 
different forms according to the kind of molecule; the virtue of 
equation (2.5) is that it gives a value for Xf{ which is valid for 
closed ring molecules as well as for simple chain and branched 
chain molecules. 

2,5. The configurational energy of the assembly is given by 


W 


2 2 2 

i i<k 


where 1^-,., 1^^ are the average interaction energies due to closest 
neighbour i—i and i—k pairs respectively. This can be written 


as 


TF 


i i<k 

~ 2 2 ^ik'^ik’ 

i i<k 


( 2 . 6 ] 

(2.7] 


where are the energies of mixing of each pair of components, 
and are defined in equation (1.13), and —x,- (all i) is the average 
potential energy of a molecule in the pure liquid i. The grand 
partition fimction of the assembly can then be written as 


Z 




( 2.8 
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where ^UV-. written for the number of way< of arrandng 

X molecules of species i (all i] on the array of sites so that there 
are zX^j closest neighbour contacts of the kind i—j tall b j): 
that is. giXj^ X,-:) is the number of distinguishable arrangements 
of the assembly of systems i (all i) a specified configura¬ 
tional energy' given by ecpiation (2.7). If all the energies of 
mix ing are zero, then the grand partition function is 


^ = Z ?(V.)expj 3 qfXiXi IT {iv,li{T]j,iT)V^<. (2.9) 


-V. 


where ) is the number of ways of arranging *V,- systems i (all 
i) on the array of sites, independently of the numbers of closest 
neighbour contacts of different kinds, since now the configura- 

C w 

tional energy depends only on the numbers of systems of different 



The essence of the determination of the grand partition func¬ 
tion is to evaluate the combinatoix' factors ) and glX^.X^j). 
For the case of a mixture of roughly spherical molecules of about 
the same size with zero energies of mixing this is trivial, it is 
simply I 2 -V, )’. U V,-'.. In aU other cases its evaluation presents 

"T i 

considerable difficulty and approximate methods must be used. 
The evaluation of is equivalent to the evaluation of an 

integral such as that given in equation (1.15). 


2. 6 . The nature of the difficulty can be seen from a considera¬ 
tion of a binary- mixture of monomer molecules. From equation 
( 2 . 3 ) it is clear that of the five quantities X,;. X.^. X. 

Xj which can be used to specify such a binary mixture only three 
are independent. Ordinarily, and as is imphed by the way in 
which equation (2.S) has been written, different variables are 
used to enumerate the complexions of the assembly and to 
determine its configurational energy. Whereas the distin- 
cniishable configurations of the assembly are specified by the 
occupation of individual sites of the array, the configurational 
enercw arises from closest neighbour (or more distant) inter- 
actions and has to be determined by the specification of closest 
neighbour (or more distant ) pairs of sites. It therefore follows 
that if a means could be found by which both the evaluation of 
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the configurational energy and the enumeration of the com¬ 
plexions of the assembly could be carried out in terms of one and 
the same kind of variable (that is, by the specification either 
of pairs of sites or of the occupation of individual sites, but 
not both), one could hope greatly to simplify the problem. The 
configurational energy can be expressed only in terms of the 
specification of pairs of sites; can, then, the distinguishable 
complexions of the assembly also be expressed in terms of these 
variables ? The usual procedure is to take N^, Nj, and Xfj as 
independent variables; the suggested procedure is to take 
Xii, Xjp and X,-,- as independent variables. This possibility has 
been examined for a binary mixture of monomer molecules, in 
particular by AHrey and Mark,t but it has been shownj that their 
treatment is based on an untenable assumption. Furthermore, 
such a specification introduces difficulties no less formidable 
than those which it was designed to overcome, so that its use 
must be abandoned. These difficulties had been previously 
recognized§ and an attempt had been made to introduce a 
correction for the mterdependence of the distribution of pairs 
of sites of different kinds. 

The nature of these difficulties can be seen from the following 
examination of the method used by Alfrey and Mark. They con¬ 
sidered a set of \N numbered pairs of closest neighbour lattice 
sites which are chosen in such a way that each site of the lattice 
belongs to one and to only one pair. If amongst the selected set 
of \N closest neighbour pairs there are X i—j pairs, then the 
numbers oii—i andj—j pairs are ^(N^—X) and X) respec¬ 

tively. The number of configurations is then calculated by con¬ 
sidering (i) the number of ways of arranging X indistinguishable 
elements of one kind, KN^—X) indistinguishable elements of a 
second kind, and ^{Nj—X) indistinguishable elements of a third 
kind amongst the available places, and (ii) by allowing for 
the fact that each i—j pair can be arranged in two distinct ways 

t Alfrey, T., and Mark, H., J. Chem. Phys. 10 (1942), 303. 

t Miller, A. R., ibid. 15 (1947), 513. 

§ Fowler, R. H., and Guggenheim, E. A., Proc. Roy. Soc. A 174 (1940), 189, 

§ 8 . 
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within its reference position. On this basis, they derive a com¬ 
binatory factor 


^ponding to this, a partition function 

J = 2 ^(X)exp(—Xu>,y/fer), 


( 2 . 11 ) 


apart from factors allowing for the internal motions of the mole¬ 
cules. 

This treatment involves the assumption that if in the selected 
group of closest neighbour pairs there are X pairs of the kind 
i—j, then in the complete assembly the total number of i~j 
pairs is exactly zX. It also implies that the total numbers of i —» 
and j—j pairs are X) and ^z{Nj—X) respectively. In 

fact, in terms of our notation, the assumption is 

zX = zXjj, ^ 

iz(A;-X) = zX,,. = iz(x;.-x,,.). (2.12) 

lz(Nj-X) = zX,.,. = ^z(i^.-X,,).J 

This assumption is essential in the treatment given by A l fr ey 
and Mark for two reasons: (i) In order that the interaction 
energy, which we have seen is equal to X^jW^p can be expressed 
in terms of the parameter X which they introduce, namely, as 
Xwp. This is the specific form in which Alfrey and Mark express 
this assumption in equation (7) of their paper, (ii) It is also 
necessary because their treatment assumes that the particular 
arrangement of particles on the lattice in which there are zX,-j, 
zX,j-, and zX^y pairs of i—i, i—j, and j—j closest neighbours 
respectively is equivalent to, and can be described precisely by, 
the aUocation of X, I(A"i-X), and U^)-X) respectively i-j, 
i—i, and j —j pairs amongst the group of |A pairs of sites selected 
in the wav specified. This is considered further in the discussion 

below (equations (2.14) and (2.15)). 

In fact, from a lattice of .V sites in which no closest neighbours 

of any site are also closest neighbours of one another, a group of 
closest neighbour pairs which are such that every lattice 
point belongs to one and to only one pair can be chosen in 
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((2— different ways.t While, at first sight, it appears 

possible that the assumption made by Alfrey and Mark might 
be vahd (or nearly so) for a sufficient number of these possible 
selections of the group of pairs to enable the method to give a 
reasonable estimate of the combinatory factor, they offer no 
justification of it. Furthermore, it is difficult to discover any 
formal justification of their assumption. When the problem is 
examined from a somewhat different point of view, the nature 
and implications of the assumption made by Alfrey and Mark 
become clearer, and it is seen to be far from justified. 

It has been noted above that A,-,-, and Xjj can be taken as 
the three independent parameters which are used to specify 
the assembly. With this specification the determination of the 
combinatory factor could then be regarded as a question of 
determining the number of ways in which zX^ indistinguishable 
elements of one type, 2A,-_^ indistinguishable elements of a .seconrl 
type, and zXjj indistinguishable elements of a third type can be 
distributed amongst izN available pairs, subject to any restric¬ 
tions inherent in the interconnexions of the assembly. The closing 
phrase indicates where the difficulty lies. A little reHection show.s 
that if from the total number of available pairs, zA",, are chosen 
for the i—i closest neighbours, then not all the remaining 
\z(N—Xii) pairs can be utilized to accommodate a j—j jjair, 
for it is inherent in the interconnexions of the pairs that a large 
number of those now remaining can accommodate only an i—j 
pair. Furthermore, the position is not eased by selecting the 
pairs for the i—j closest neighbour pairs second; for unles-s severe 
restrictions on our freedom of choice of the required number of 
pairs are accepted, w'e would be left finally with a large number 
of pairs which could not accommodate a J— j closest neighbour 
pair and at the same time an insufficient number of possible 
pairs to accommodate all the j—j closest neighbours. In pas.sing 
it should also be noted that the precise specification of tliese 
restrictions, if not impossible, is certainly neither simple nor 
straightforward. It appears likely that it is in this restriction 
on the free choice of elements to accommodate the closest 

t Miller, A. R., Proc. Cambridge Phil. Soc. 38 (1942), 109, equation (16). 
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neighbour pairs of different kinds that the implications of the 
assumption made by Alfrey and Mark are to be found. We now 
proceed to show that this is, in fact, the case. 

Stirling’s theorem for the asymptotic behaviour of a factorial 
ensures that for large values of r and s and a constant a we have 


a log 


(r+s)l 

r!s! 


1 Wr+fi’. 

^(ar)!(a5)! 


(2.13) 


If the assumption made by Alfrey and 
can be justified, it then follows that 


!ark, equation (2.12), 


(m 




(izN)! 


The right-hand member of equation (2.14) reduces to 

(M)'- 


(2.14) 


(zx,,y.(zx,jr. (zx^,)r 


(2.15) 


which is precisely the number of ways of arranging zX^ elements 
of one kind, zXf^ elements of a second kind, and zZ,-,- elements 
of a third kind amongst ^zN positions vnthout any restriction 
whatsoever. Thus the assumption made by Alfrey and Mark 
amounts to neglecting the restrictions on the allocation of the 
closest neighbour pairs of different kinds which are inherent in 
the arrangements of two or more species of molecules on a regular 
lattice. Thus it gives a value for the combinatory factor which 

is too large and far from correct. 

When the molecules of one or more species of the mixture con¬ 
sist of chain molecules, the number of closest neighbour pairs 
occupied by adjacent sub-molecules must also be specified. 
These are the quantities defined by equation (2.5). In the 
general case the suggested procedure would be to take quantities 
of which Xii, Z„, and are representative as independent 
variables. Here there is also an interdependence between the 
ar pairs of closest neighbour sites which are occupied by adjacent 
segments of the polymer molecule. In short, the x^ pairs have 
to be selected according to the sites which can accommodate 
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one of the polymer molecules; and there is no simple way by 
which this constraint can be taken into accoimt. It is a constraint 
of precisely the same kind as that implicit in the specification ol 
the numbers of closest neighbour contacts of different kinds, 
which has been dealt with in detail in the preceding paragraphs. 
Furthermore, this difficulty persists for mixtiu^s for which the 
energies of mixing are zero. We conclude that no simple distri¬ 
bution of pairs of sites of different kinds amongst the available 
pairs of sites on the array can provide a means of evaluating 
the combinatory factor in any of the cases of interest. Recourse 
must therefore be made to indirect methods of evaluating the 
relevant partition functions. We now proceed to examine these 
methods. 



Ill 

EVALUATION OF THE COMBINATORY FACTOR 


3.1. To evaluate the partition function for a co-operative 
assembly we have seen that it is necessary to proceed by some 
indirect method. U the equilibrium properties of the assembly 
can be determined by an indirect method, then it will clearly be 
possible to evaluate the partition function, and, in particular, 
the combinatory factor by a process which is just the reverse 
of that normally used to determine the equilibrium properties 
from the partition function. The purpose of this chapter is to 
consider some of the methods which have been used to determine 


the &ee energy and the equilibrium properties of a co-operative 
assembly, and from them to determine the combinatory factor 
which occurs in the partition function for such an assembly. 
The treatment is confined to the scope of this monograph and 
is not intended as an exhaustive examination of the statistical 
mechanics of co-operative assemblies in general. Attention is 
focused on those methods which have proved most useful in the 

theory of solutions of high polymers. 

The free energy of a co-operative assembly can be expanded 
formally as a power seriesf in l/kT, in which the coefficients of 
the successive terms are the cumulants of the unweighted 
averao^e values of the moments of the configurational energy. 
In effect, this gives a power series expansion for the quantity X 
defined in equation (1.30), and the moments have been evaluated 
as far as the third. In this way it can be shown that the approxi¬ 
mation which had been used by GorskyJ and by Bragg and 
Williams§ in considering binaiy alloys was equivalent to neglect¬ 
ing all but the first term in the expansion, and that the method 
introduced by Bethe.H in which local configurations are con¬ 
sidered in detaU, is vahd up to the third moment of the con- 


Kirkwood, J. G., J. Chem, Phys. 6 (1938), 70. 

+ Gorsky', W., Zeit.f. Physikj 50 (1928), 64. .oojv can 'Ki/i 

s W T. find Williams. E. J.. Proc. Roy. Soc. A 145 (1934), 669; ibid 



Ch.in,§3.1 EVALUATIOX OF THE COMBIXATORY FACTOR 29 

figurational energy, f It has also been shownj that to examine 
the properties of an extensive class of co-operative assemblies 
the method of chemical kinetics or the condition of quasi¬ 
chemical equilihrium§ is equivalent to Bethe's method in its 
first approximation. Since an unambiguous formula for the com¬ 
binatory factor, for a mixture containing pohmier molecules 
with zero energy of mixin g, was first obtained by using Bethe’s 
approximation, we wiQ now examine this approximation and 
its use in this connexion. 

Bethe’s method 

3.2. In the Bragg-WiUiams approximation the evaluation 
of the integral in equation (1.15), or X in equation (1.30), or of 
g{Ni, Xfj) in equation (2.8) is avoided by assuming that the force 
tending to produce order at a given point is determined uniquely 
by the average state of order throughout the assembly. Actually 
the order at a particular site of the array depends on the con¬ 
figuration of the particles in its immediate neighbourhood and 
only the average local configuration is determined by the long- 
range order of the array as a whole. There are fluctuations about 
this average configuration and corresponding fluctuations in 
the configurational energy. The neglect of these fluctuations 
simplified the problem by making a less precise specification of 

energj" adequate. To take 
the fluctuations into account, it is necessary to determine X or 
g{Ni, Xfj), as will have been evident from the considerations of 
§§ 1.15 and 2.5. 

An approximate value of the partition function taking account 
of the local fluctuations can be obtained in the following way.H 
It is simplest to outline the argument for a binary quasi-crystal¬ 
line mixture of simple molecules. The occupation of a small 
group of sites is considered in detail. Suitable groups are one 

t Bethe, H. A., and Kirkwood, J. G., J. Chem. Phys. 7 (1939), 578. 

t Rn^brooke, G. S., Proc. Boy. Soc. A 166 (1938), 296; Guggenheim, 

E. A., ibid. A 169 (1938), 134; Fowler, R. H., and Guggenheim, E. A., ibid. A 
174 (1940). 189. 

§ See § 5.3. 

i Bethe, H. A., Proc. Boy. Soc. A 150 (1935), 552. 


states of different configurational 
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site, such as site 0 in Fig. 1, a pair of closest neighbour sites, such 
as sites 0 and 1 in Fig. 1, a central site and its first shell of z 
closest neighbour sites, such as sites 0, 1, 2, 3, 4 in Fig. 1, and so 
forth. In the first case the assembly consists of one site (the 
central site) and a residue which forms an external field of force 
so far as the central site is concerned; in the second it consists 

• • • • • 



Fig. 1. Group of sites in applying Bethe’s method 
(drawn for the two-dimensional case). 

of a single pair of sites and a residue which provides an external 
field of force; in the next case the assembly consists of a central 
site and its closest neighbours {z -\-1 sites in aU) in a field of force 
provided by the remaining N—z—\ sites; and likewise for larger 
groups. In considering a group of 2 +1 sites consisting of a 
central site and its first shell of closest neighbours we impose 
the further restrictiont that no closest neighbours of a given site 
are also closest neighbours of one another. The effect of this 
assumption has been examined for regular solutions of simple 
moleculesj and for solutions of polymer molecules.§ Oose- 
packing of the molecules is found not to introduce any significant 
difference in the results. Accordingly, because the simpler 
algebra makes the nature of the processes clearer, this restnction 
wiU be retained in this exposition of the method. The mter- 

t Cf. Bethe, H. A., Proc. Boy. Soc. A 150 (1935). 552. 

% Buahbrooke, G. S., ibid. A 166 (1938), 296. 
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ca». in. § 3.2 - 

actions of the particles which occupy the selected group of sites 
are into account explicitly and equiUbriuin conditions are 
obtained by the requirement that both central and first shell 
sites must be average sites. Factors, e, are used to specify the 
occupation of these sites in the foUoiidng way: factors c^, are 
introduced for each molecule of type i which occupies the central 
or a first shell site respectively; factors e,-, ^ are introduced for 
each molecule of type j which occupies the central or a first shell 
site respectively. and are simply related to the corre¬ 
sponding absolute activities; ej and ^ include in addition factors 
to take account of their interactions with particles in the residue 
of the assembly. The grand partition function is then given by 




(3.1) 


where tj is written for exp{—WiJzkT) with (as before) the 
energy of mixing. 

3.3. Some light is thrown on the nature of Bethe’s method by 
a variant of it proposed by Guggenheim.f The grand partition 
fdnction for a binary mixture is (cf. § 2.5) 

Z = 2 (3.2) 

where the partition function of a molecule of type % for all degrees 
of freedom and all forms of energy other than the interaction 
energy has been denoted by Thus 


Z 


2 m. ifj, 

Ni,Nj 


(3.3) 


where has been written for exp(—i<7^y/zfe!r) and for 

respectively. 

In equation (3.3) the series for Z can be replaced by its maxi¬ 
mum term without introducing other than negligible errors in 
deriving averages. If the values of N^, and in the maxi¬ 
mum term be denoted by asterisks, then for the mean values we 

»i = Nt, N, = Nf, X, = Z*, (3.4) 


and 


Z = Z* 


j(Zf. Nf, 


(3.5) 


Now define geometrical means in the following way. Introduce 


Chigganheim, 



32 EVALUATION OF THE COMBINATORY FACTOR Ch. Ill, §3.3 

a quantity ^ which, is the geometrical mean of the factors of the 
type contributed by each site of the array to Z*; a quantity tj 
which is the geometrical mean of the factors of the type rjy 
contributed to Z* by each pair of closest neighbour sites; a 
quantity y which is the geometrical mean contribution per site 
to g{Nf,Nf,X*), and quantities typified by 17 ^ which is the 
geometrical mean of the factors contributed to Z* by each 
pair of closest neighbour sites of which one (at least) is specified 
as being occupied by a molecule of type i. Thus 


= If 




izN 


Vii 


(3.6) 


yN=g{NlNf,X*). 


Now consider a central site and its z closest neighbours, and write 
down the terms contributed to Z* by these z+1 sites explicitly 
and use the averages given by equations (3.6) for the contribu¬ 
tions to Z* made by the N-z- 1 sites forming the residue of the 

assembly. This gives 

y^yN-Z-l^N-Z-ly^izN-^^ ( 3 . 7 ) 

where the factors of the type Tjf-i are introduced to aUow for the 
interactions between particles which occupy first shell sites and 
the outer sites. The exponent of the factors of the typei^j is neces¬ 
sarily z -1 since of the z closest neighbour interactions enjoyed 
by each particle which occupies a first sheU site one, that with 
the central site, is introduced expUcitly through the factors 17 ^^. 


write 


we obtain 


1 /= 


(3.8) 


£—1 




V VW A* A 


(3.9) 

Since we are concerned with ratios of tie terms in equation (3.9). 
which give the relative probabilities of different coua^ratiom 
it. IS dear that the common factor T 
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omitted without possibility of error. This factor merely gives an 
average contribution from all sites other than the z-{-l sites in 
the selected group. Thus correct results will be obtained by 
writing 




(3.10) 


for the partition function. This is precisely the expression which 
is obtainedf by constructing the grand partition function accord- 
ing to the method which was proposed by Bethe. The identity 
of equations (3.1) and (3.10) indicates the nature of Bethe’s 
approximation and the physical significance of the e factors. 
It requires no elaboration to establish that for a binary mixture 
it is only necessary to introduce e factors for one of the two 
types of molecules, for the significant result will be the value of 
the ratio ejcj and not the individual values of and e^. Generally, 
however, in the sequel we shall be concerned with the ratio of 
the €i, €j to the ‘e-factor’ appropriate to a vacant site. As will 
be seen below, admitting the possibility of imoccupied sites 
provides a check on the internal consistency of the method. J 
The way is now clear to present the discussion of mixtures con¬ 
taining dimer and trimer molecules§ which led to the formula¬ 
tion of an expression for the combinatory factor appropriate to 
mixtures containing polymer molecules containing any number 
of segments. 

Mixtures containing dimer molecules 

3.4. A mixture of monomer and dimer molecules which 
can be arranged on a quasi-crystalline array of sites is considered. 

There are N sites in all, and 

# 

N = N^-^2N,+N„ 

where is the number of sites left vacant when all the monomer 
and dimer molecules have been arranged on the array of sites. 
The occupation of a group of sites consisting of a central site 
and its z closest neighbours (none of which is a closest neighbour 
of any other) will be examined in detail. The fractions of the 


t Equation (3,1) above. 

t Miller, A. R., Proc, Cambridge Phil. Soc. 42 (1946), 303. 
§ Miller, A. R., ibid. 38 (1942), 109; ibid. 39 (1943), 54, 
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total number of sites which are occupied by monomer and dimer 
molecules respectively wUl be denoted by and d,, so that 

In applying Bethe's method the parameters will be defined in 
the following way. For a monomer, and wiU be introduced 
according as it occupies the central or a first shell site. For a 
dimer which occupies the central site and a first shell site a para¬ 
meter e, will be introduced, and for a dimer which occupies a 
first shell site and an outer site a parameter e'j. It is worth while 
to consider this case in somewhat greater detail than will be 
necessary in considering other cases subsequently. Let the 
number of monomer and dimer molecules occupying the central 
site be denoted by ^nd respectively; each of these can 
assume the values zero and unity, but if either is umty the other 
is necessarily zero. Further, let the numbers of monomer and 
dimer molecules which occupy first shell sites be denoted by 
and n.-> ^ respectively. The weight factor, that is, the number of 
distinct ways in which a particular configuration specified by a 

given set of values of n^ o, n. 2 , 0 , ” 2,1 ^ 

denoted by g. All the ways in which the group of sites consisting 

of the central site and its closest neighbours can be occupied are 

indicated in Table 1, which for the sake of clarity is given for 

the case in which 2 = 4. Appropriate tables for other values of z 

can be constructed in a similar way. For instance, when is 

zero and n 2 ,o is unity, the sum of and n .21 cannot exceed z-l; 

can be given the set of values 0, 1 , 2 ,..., z—1; if has the 

value r (0 < r < z-l), then can assume the set of values 

A 1 9 z-r—1 The number of distinct arrangements m 

which the set of particular values 


n 


0 . 




1 


10 — ■• 2,0 ' 

where 0 ^ ^ 2 —1, 0 ^ ^ 


72 


1,1 




n 


2,1 


2 ’ 


r^— 1 , can arise is given by 




The corresponding weight factors and range of allowable values 
of and can be determined in a similar way for the othei 

values of Hi,o und n.y_o- 
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Table 1 


Occupation of the group of z -\-1 sites 


(a) nj , zero, unity. 



(t>) fhji, zero, zero and (c) unity, n*zero. 



The typical term in the grand partition function constructec 
according to Bethe’s method is 


Here t) has been written for esp{—w^^zkT), where 


is th( 


energy of mixing of monomer and dimer molecules, and th< 


energy zero is chosen so that the total interaction energy oj 
the pure species is zero. This term gives the relative probability 
of that mode of occupation of the group of sites which is specified 
by the given set of values of and Wa j. It therefore 

follows that if the partition function is constructed according 
to Bethe’s method we getf 



e2)^+€i{l+ci+i]€2)^-l-2e2(l+Tjei+e2)*-i. 



3.5. The frequencies of occupation of the central site by mono¬ 
mer and dimer molecules respectively are given by c^Slog Zjdey 
and € 2 Slog Zjde.^. Since the central site is chosen at random its 
frequency of occupation by different molecules must be the 


t The corresponding result for the two-dimensional problem of an adsorbed 
mMolayer was given by Roberts, J. K., and Miller, A. R., Proc. Cambridgt 
Phil. Soc. 35 (1939), 293, and by Chang, T. S., Proc. Boy. Soc. A 169 (1939), 
W2, and for a three-dimensional array of sites, aU of which are occupied, by 
l^er, A. R., Proc. CambridgePhU. Soc. 38 {1 942), 109. It should be noted that 
the treatment given by Chang in § 3 of his paper is open to objection, as is also 
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le as that of any other site, that is, it is an average site, 


ence 






e. 




(l+ei+c;)- ’ 

2^2(1+'>?cl+ 4 )^ 


(3.12) 


(3.13) 


len the central site is unoccupied, the average numhers of 
;t shell sites which are occupied hy monomer and dimer mole- 
es respectively are given by e^dlog ZJde'j^ and e'gSlogZo/Seg, 
lere has been written for the first term in the right-hand 
jmber of equation (3.11). When the central site is occupied hy 
nonomer, the average numhers of first shell sites which are 
Bupied by monomer and dimer molecules respectively are 
ren by e^dlogZjJdei and e^dlogZJde^, where Z^ has h^n 
itten for the second term in the right-hand member of equation 
.11). Likewise, when the central site is occupied by a dimer, 
e average numbers of first shell sites which are occupied 
r monomer and dimer molecules respectively are given by 
0 log ZJd€\ and l+e^8 log ZJde'i, where Z^ has been written for 
le third term of the right-hand member of equation (3.11). 
nee the first sheU sites are average sites their overall average 

nonomer and dimer molecules respectively must 

Thus 


:cupation by 
> z6i and 


td, = (1-0X 


ze, = (1 




,8logZ 



hich the energy of mixin g 
reduce to 


01 


1 - 01-0 


(3.16) 
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In this case also, equations (3.14) and (3.15) give 

z 6 i{l-\-ei-{-e 2 ) = (z—(3.1! 
2^2(l'f'^l"t"^2) = (3.1i 

The average number of vacant first shell sites is given by 



whence z(l—0^— 02 )(l+€i+e 2 ) = z— 02 - (3.2( 


By adding equations (3.18) and (3.19) and subtracting the sm 
firom 2 ( 1 + 61 + 62 ) equation (3.20) is recovered as is necessar 
Substituting the value of I+ 61+62 given by equation (3.20) i 
equation (3.17) we obtain 



Equations (3.16) and (3.21) give the essential results for 01 
immediate purpose. 


Partial differential equations for the 

combinatory factor 

% 

3.7. It has been seen (§ 3.3) that the parameters 6^ which ai 
introduced in Bethe’s method have to be identified with 
which is the product of the absolute activity and the partitic 
function for all degrees of fieedom and all forms of ener^ 
other than the interaction energy. The relation between the 
and the combinatory factor can be found by considering tl 
grand partition function. For a mixture of any number of speci 
it can be written in the form 




EVALUATIOX OF THE COMBINATORY FACTOR Ch. Ill, § 3.7 


jilibriiim values A X^j of the A,- and the are given by the 
s of equations 

r' 

VlncTr7/Y. Y 


cX. 


^ log g{Xf, X,,.)+log = 0 (aU i), 


(3.23) 




log giXi. Xfj )-f 2 log = 0 (all i, j). 


(3.24) 


order to obtain an expression for the corabinatoix- factor it is 
iumed that the eqiuhbrium conditions determined by Bethe's 
?thod are precisely equivalent to those determined by the sets 
equations (3.23) and (3.24). That is, that the 0,- (determined 
equations such as (3.16) and (3.21)), u hich give the average 
ictions of the total number of sites which are occupied in a 
ecified manner, and the iV) determined by equations (3.23) 
d (3.24) give identical specifications of the equilibrium state. 


ms 


Xe, = r,iY,, 


lere r,- is the number of sites which are occupied by each mole- 
Je i. 

For the case in wliich the energies of mixing are zero there is 
nply a set of equations corresponding to (3.23). Identifying 
e (i with the e,- of Bethe's method and writing the equations 
terms of the 6^ rather than in terms of the X^, these give, in 
neral, for a mixture of two species of molecules, 


r,. e 


X ee.. 


log?(.^;•,i^■;•,.^;)^-log€,• = o. 


(3.25) 


^ ^ log V)^log «;• 


0 . 


1 order that these partial differential equations shall be in- 
■grable it is necessary- that 


9 ' 


log e,- 


. — log e,-. 
’c0.- ' ' 


(3.26) 


nd the integrabihtv of the partial dififerential equations pro- 
ides a check on the internal consistency and the vahdity of 


A ^ Vh ^ f\ 


• 1 * ^ 
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Combinatory factor for a mixture containing 

dimer molecules 

3.8. From the values given for Cj and eg by equations (3.16 
and (3.21) it can be verified that 






dd^ 


log eg. 


Thus the partial differential equations 


1 8 

N801 


log9(^i,-A^2.-^t.)+log^i-log(l-0i-02) = 0, 


N 88. 


logg{Ni, N^, i^J+log 02 +log(z- 02)-2 log(l -^i-^g)- 

— 2 log 2 


are integrable; and integrating them the combinatory facto 
g(Ni,N^,NJ can be evaluated. The constant of integration i 
determined by the fact that when is zero, g{Ni, N^, N^) mus 
reduce to the number of ways of arranging indistinguishabl 
particles on sites leaving N„ sites vacant, that is, by th 

fact that g(N^,0,N,) = (N,+N,)\IN,\N,l 

Integrating the above partial differential equations we obtain 


logy(i'^i,-2V^2,-Z'4) 


01 log 01 —| 02 log 02 +^ 2 —^ 2 )log (z—^ 2 ) 


— (1 —01—02)log(l —01—02) + 02log2—^Zl0g2. ( 3 . 2 ': 

If 01 and 02 are replaced by their values in terms of .^i, N^, an 
K we obtain finally 


g(Ni, N^) 


{N,+2iz-l)NJz+N,}[VHIi 


{N^+2N^+N,) 


K}\yHN^+2N,+N,)\lzY 

! N,\N,\NJ W 


(3.21 


where a is the symmetry number for the dimer molecules and: 
equal to 2 if the ends of the molecule are indistinguishable, an 
to unity otherwise. 


Mixtures containing trimer molecules 

3.9. A mixture of monomer and trimer molecules whic 
can be arranged on a quasi-crystaUine array of sites is considere* 
There are N sites in all, and 
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where is the number of sites left vacant when all the monomer 
and trimer molecules have been arranged on the array of sites. 
As in the case of the mixtm-e containing dimer molecules, the 
occupation of a group of sites consisting of a central site and its 
first shell of s closest neighbours is examined in detail. The 
fractions of the total number of sites which are occupied by 
monomer and trimer molecules respectively will be denoted by 
01 and 03 - The Bethe e-parameters will be defined in the following 

lor a trimer which has its centre sub-molecule on the 
site of the group a parameter €3 will be introduced; for 
trimer which has an end sub-molecule on the central site of 
le group a parameter €3 will be used, and for a trimer which 
IS an end sub-molecule on a first shell site and its other sub- 
olec\iles on outer sites a parameter €3 will be used. Parameters 
, cj will be introduced for the monomer molecules according as 
ley occupy the central or a first shell site. The detailed occupa- 
on of the selected group of sites can be considered in exactly 
le same way as was done in § 3.4 for mixtures containing dimer 
olecules. When the energy of mixing is zero, Bethe’s approxi- 
ate grand partition function is 



+ Jz€3(l-i-€i+e;)»-2 (3.29) 

= •^0 + + Z ^ Z ^, 

which serves to define the quantities Z^, Zj, Zg, and Z 3 which 
are introduced for convenience of notation. Equation (3.29) 
is the grand partition function for the case in which the trimer 
molecules are completely inflexible. The case of perfectly 
flexible molecules can also be considered without difficulty, for 
the only difference is that the numerical factor ^2 in the quantity 
Z 3 must be replaced by ^ 2 ( 2 —1). This makes no significant 
difference in the analysis. 


3.10. The frequencies of occupation of the central site by 
monomer and trimer molecules respectively are given by 
€i0logZ/0ei and Z jee^-^-e'^dloQ Z Since the average 
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mode of occupation of the central site must be the same as that 
of any other site it follows that 


01 


i-e.-9 


( 3 . 30 ) 


and 


0. 


i-d,-0 


1 ^3 

ze3(l+e^-f~g3) + l^^3 


( 3 . 31 ) 


In those configurations in w 
trimer molecule it must be ( 
as many configurations as i 
of a trimer. Applying this i 
of the group it follows that 


2c3^1ogZ-4^1ogZ, 


Be 


Be. 


whence 


^3 — 


( 3 . 32 ) 


Now consider the occupation of the first shell sites by monomer 
and trimer molecules. Considering separately the occupation of 
the first shell sites according to the mode of occupation of the 
central site, in the same way as was done in the case of a mixture 
containing dimer molecules, it follows that the frequencies of 
occupation of the first shell sites by monomer and trimer mole¬ 
cules respectively are given by 


Be 


Be- 


and 


Be^ 

Bei 


+¥>[ 2 + 4 ^-^] 


These expressions have the values 


( 1 - 0 ,-03) 


Z€ 


l+el+e; 




Z€ 


l+ei+ej 


1 + €i + €8 


-Llfl 
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ind 




Z€ 


l+el+e' 


» + ^l 




1 + 


(z 


1+4+ 



+ 


+i^3f 


2 + 


(z-2)6; \ 

i+4+4/‘ 


(3.34) 


a 

LJ 


lince the first shell sites are average sites, expressions ( 3 . 33 ) and 
3.34) must be equal to z0j^ and z0^ respectively. This leads to the 

quations z 0 ,(l+,;+e;) = z.;+|(l+ 4)03 ( 3 . 35 ) 


and 


2 ^ 3 ( 1 +4+4) = 2€'i—| e'i 0,. 


(3.36) 


By analogy with expressions (3.33) and (3.34) the number of 
vacant first shell sites must be given formally by 




1+4+4 


i+^l 


1+4+4 


+ t^3T";— r-,—i + 


1+4+4 




z —2 


which leads to the equation 


1+4+4’ 


z(l — 01 —03)(l + fl + ^3) 


2—503- 


(3.37) 


If we add equations (3.35) and (3.36) and subtract the sum from 
s(l+ei+€ 3 ) we recover equation (3.37) as is necessary. Equa¬ 
tions (3.31), (3.32), and (3.37) can be solved for €3 and we obtain 


2 03 ( 2 - 103 )^ 

23(1-0^-03): 


(3.38) 


The partial differential equations for N^, N^) are therefore 


1 d 


Nd0i 


Iog 0 (i 4 ,i\r 3 ,.vj+iog 0 i—iog(i— 01 — 03 ) = 0 , 


3 d 


Nd0. 


log 0 (A\, iVg, iv;,) +log 03+2 log(z— 103 )—3 log( 1 — 01 — 03 )+ 


+log 


3z3 


0. (3.39) 


From the values given for Cj and €3 by equations (3.30) and 
(3.38) it can be verified that 


3^1ogCi = 3(1-01-03) 


-1 




log €3. 
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Thus the partial differential equations (3.39) are integi 
Integrating them we obtain 

ilog N 3 , IQ 

= -01 log 0 i-i 03 log 03 -(l -^i-^3)log(l -^ 1 -^ 3 )+ 


+ K 2 -i^ 3 )l 0 g{z-i^ 3 )+i ®3 log 


3z3 


^zlogz, (3.40) 


where the constant of integration has been determined by the 
fact that when (that is, N^) is zero g{N^, N^, N^) must reduce to 

K 01 and 6 ^ are replaced by their values in 

terms of .^ 1 , N^, and we obtain finally 

.r. r{^i+(3z-4)i\r3/z+i\^}niMi'^l+3i^3+^»)!/^^* 

0 (^ 1 ,iVg,ivj - (iVi+3JV3+i\rj! J N,\N,\NJ W ’ 

(3.41) 

where a is the symmetry number for the trimer molecules. 
Equation (3.41) applies to completely inflexible molecules, that 
is, each trimer molecule is constrained to occupy a linear group 
of sites. 

3.11. An estimate of the effect which this assumption has on 
the result can be obtained by considering the other extreme 
case of perfectly flexible molecules. It has been mentioned (§3.9) 
that the only effect which this has is to alter, by a multiplicative 
factor s— 1 , the numerical factor in one of the four terms of the 
grand partition function (equation (3.29)). It is easy to see that 
this will merely multiply the term which involves £3 in equations 
(3.31) and (3.32) by a factor z— 1 . Corresponding to this the 
value of ^3 is altered, in an equaUy trivial manner, to give 


2 eJz-USo)^ 
32»(z-1) (1-01-03)®’ 


(3.38') 


and the combinatory factor for perfectly flexible molecules is 
g{N„N„N„) 

{N^z-^)NJz+N„}\-\i- rz(z-l)i^* 


{N,-\-ZN^+N,)l 


i± 3 ^±^!r; 

N^IN^INJ L 


(3.41') 
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Thus the perfect flexibility of the trimer molecules introduces, 

into the combinatory factor, a factor 2 — 1 for each trimer mole¬ 
cule. 


3.12. To consider solutions of longer chain molecules by 
exa minin g the detailed occupation of a central site and its z 
closest neighbours is excessively complicated. However, pro¬ 
vided one is concerned with systems with zero energy of miTing 
it is sufficientf to examine the detailed occupation of a group of 
sites just large enough to accommodate one molecule of any 
species in the mixture. Thus provided one considers only rigid 
chain moleciiles it is sufficient to examine the detailed occupa¬ 
tion of a collinear group of sites equal in number to the number 
of sub-molecules in the chain molecule. In this way Guggen¬ 
heim has considered mixtures containing rigid tetramers. His 
analysis leads to results which are consistent with equation 
(3.48) below. 


The free energy of mixin g 

3.13. The contribution to the free energy of a mixture which 
arises from the configurational term in the partition function is 
kT times the natural logarithm of the appropriate combinatory 
factor. Thus for a mix ture of two kinds of roughly spherical 
molecrdes of about the same size the configurational contribution 
to the free energy of a mixture containing of type i and of 

type j is ^ „ 

-kTNM^-kTN,^^. (3.42) 


For a mixture which contains dimer and monomer molecules 
the configurational contribution to the free energy isj from 


equation (3.28) 
k T{^zNi -f (3—1 )-^2}log 


N,+2{z-l)NJz 

N^+2N, 


kTNilog 


N,+2N, 


Expression (3.43) reduces to zero for a sample of the pure liquid 


t Guggenheim, E. A., Trans. Faraday Soc. 41 (1945), 107. 

1 We the caae obtained by putting zero in equation (3.2$). 



(f monomer molecules, when is set equal to zero. On the 
ither hand, the pure liquid of dimer molecules has a combinatory 
erm in its free energy, as can be seen by setting equal to zero 
n (3.43). This is due to the distinguishable ways in which the 
limer molecules can pack amongst themselves. For the con- 
igurational term in the free energy of the pure hquid of dimers 

TO get 

_ ___ Z —^ 1 « «« ^ ^ 2z J A A \ 


kT{z—l)N2log 


-^kTN^log 


(3.44) 


Che configurational term in the free energy of mixi ng of a mix- 
lure of dimer and monomer molecules is therefore given by 
lubtracting (3.44) from (3.43). This gives 


(A^): 

kT 


{izi\ri+(z-i)J\rjiog 


N^+2(z-1)NJz 


2^1 log 


N^+2N, 




1 )i\4 log 


1 ’ 


(3.45) 

rrespective of the value of the symmetry number a. 

The configurational term in the free energy of a mixture con- 
aining trimer and monomer molecules can be found in precisely 
he same way from equation (3.41) or (3.41'). Irrespective of 
whether the trimer molecules are flexible or inflexible, it can be 
hown that the free energy of mixing is given by 


(A4)3.i 

kT 


{i2i\ri+(3z-4)i\r3}log 


JVi+(3z—4)iV'3/2 






or all values of the symmetry number. 


(3.46) 


Mixtures containing polymer molecules 

3.14. Comparison of expressions like (3.28) and (3.41) for 
-he combinatory factor, or of (3.45) and (3.46) for the free energy 
)f mixing, of mixtures containing dimer and trimer molecules 
espectively led to the construction of the combinatory factor 



fipee 


mixture contained polymer molecules of any number of sub¬ 
molecules. The generalization is, in fact, straightforward and 
merely involves writing down the simplest general formula to 
which those which had been evaluated for mixtures containing 
dimer and trimer molecules conform. Inspection of these formu¬ 
lae shows that the numerical factors involved are in every case 
the particular values of the parameters q (defined in § 2 . 2 ) for 
the molecule in question. Thus introducing q’^, 3 ' as the values of 

q for chain molecules of two and three segments respectively, 
equations (3.45) and (3.46) can be rewritten as 


kT 




N,log 


K 


N^-\-2K 


(AA) 


3.1 


WiN^^ogq^, (3.45') 




kT 


- N, log- 


■^3 log 


K 




Izjs iyTg log 33 . (3.46') 


Once the occurrence of the values of the parameters q had been 
noticed, the free energy of mixing for a mixture which contained 
polymer molecules each consisting of r sub-molecules was given 


as 


kT 




^rlog 


N. 


N,+rN, 


WrK^ogq",. (3.47) 


responding formula 
containing rigid, rod 


K, N,) 


(N,+q',N,+N,)\y^ (N,-\-rK^NX 


{K+rK+K)\ 





(3.48) 


MiOer, A. R., Proc. Cambridge PhU. Soc. 39 (1943), 54, § 3. 



It has been seen that the perfect flexibility of trimer molecules 
introduces an additional factor z—1 per trimer molecule; for 
longer chains the flexibility wiU introduce a more complicated 
factor. It will in all cases introduce a factor per chain molecule 
R-hich will depend on the number of sub-molecules of which the 
polymer is built up, the geometry of the arra\" of sites (through 
the average coordination number) and the degree of flexibihty 
enjoyed by each polymer molecule. The average contribution 
per polymer molecule to the combinatory" factor arising, through 
the three factors mentioned, from the flexibihty of the pohuner 
molecules can be specified by a quantity Then, for a mixture 
containing r-mer molecules of any degree of flexibihty, the 
combinatory factor can be written formahy as 



(3.49) 


3.15, In evaluating the free energy of mixing it is necessary 
CO consider not only ^(.Ai,JN^) but also its value when is set 
equal to zero, which can be denoted by gr(0, .^^), and to determine 



9 ( 0 , ' 


(3.50) 


[t is clear from the form of equation (3.49) that this quantity is 
lecessarily independent of the value of Thus equation (3.47) 
'or the free energy of mixing is vahd whatever the degree of 
iexibihty of the polyuner molecules. It wih be seen in the next 
chapter that in determining the vapour-pressure equations, 
md thence aU the thermodynamical properties of the mixture, 
:he only significant quantity which does not involve the ratio 
;3.50) is the differential coefficient clogg{A\,N^}ldXj^; and it also 
nust be independent off,. Thus the quantities 'with which we 
ire concerned in determining the free energy^ of mixing and the 
rapoTir-pressure equations are necessarily independent of anv 
actor which depends on the flexibility of the polyuner molecules, 
oro'dded only that they are as flexible in solution as in the pure 
species. Since the flexibility necessarily depends on the steric 



properties of the molecules this condition will probably be ful¬ 
filled for most solvents. Thus it is worth emphasizing that since 
the precise form of is immaterial so far as the interesting 

thermod^mamical quantities are concerned, it is unnecessary to 
evaluate it. 

Thus all the quantities necessary to determine the thenno- 

d\uiamical properties of solutions containing polymer molecules 

of any number of segments have now been evaluated. Before 

proceeding to their evaluation we may first consider other 

methods which have been used to derive the combinatorc 

factor, and secondly, present a stricter derivation of the general 
equation (3.49). 

Direct enumeration of distinct configurations 

3.16. It is clear that in principle the combinator}- factor 
could be obtained bv the direct enumeration of the distinguish- 
able configurations of the molecules of the mixture on a quasi- 

crvstalline arrav of sites. To do this it would be sufficient to 

^ % 

consider the process by which polymer molecules are added in 
succession to the empty array of sites and to coirnt the number 
of alternative configurations available to each successive mole- 

O 

cule. This is a method which has been usedj to allow for the 
differences in the size and shape of the two kinds of molecules 
in a mixture of dimer and monomer molecules. In this way, 
Fowler and Rushbrooke obtained formulae in the limiting cases 
of highly concentrated and extremely dilute solutions. Their 

o V ^ 

formulae are asymptotically equivalent to equation (3.28) 
derived by Bethe’s method; a fuller comparison of the resxilts 
has been given elsewhere.! The general case of a mixture con¬ 
taining polymer molecules of any number of sub-molecules has 
been treated in a similar way by Flory§ and by Huggins,|j a 
brief account of whose work will now be given. 

t Fowler, R. H., and Rushbrooke, G. S., Trans. Faraday Soc. 33 (1937), 
1272. 

X Miller, A. R., Proc. Cambridge Phil. Soc. 38 (1942), 109, § 3. 

§ Florv, P. J., J. Chem. Phys. 9 (1941), 660; ibid. 10 (1942), 51. 
li Huggins, M. L., ibid. 9 (1941), 440; J. Phys. Chem. 46 (1942), lol; 

Ann. S.Y. Acad. Sci. 43 (1942), 1. 
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3.17. Floiy uses a model essentially the same as the quasi¬ 
crystalline array of sites (or of cells) which has been described 
in §§ 1.9, 2.1, and 2.2. He considers a mixture containing 
monomer and r-mer molecules which can be arranged on an 
array conta inin g N (= N^-\-rN^) sites in all. Starting with aU 
the sites empty, the successive addition of polymer molecules is 
examined. After n polymer molecules have been accommodated 
on the array, there are N—nr sites which may be occupied by 
the terminal segment of the next polymer molecule to be accom¬ 
modated on the array of sites. For a random mixture the average 
number of vacant sites which are closest neighbours to the site 
which is occupied by the terminal segment of the (n+ l)th poly¬ 
mer molecule is z{N — m)IN, and this is the number of sites 
available for the second segment. For the third segment, the 
average number of available vacant sites will be {z — l)(iV'— rn)jN, 
since one of the sites which are closest neighbours to that oc¬ 
cupied by the second segment is already occupied by the terminal 
segment of the polymer. Flory assumes that, on the average, 
the number of vacant sites available for every segment beyond 
the third is also equal to {z —l){Ar— m)jN, This assumption 
does not allow for the fact that, for^these segments (for explicit¬ 
ness, say, segment i), one of the closest neighbour sites of the 

{i— 1 )th segment may already be occupied by one of the segments 

0, 1 , 2,..., i—4. 

Thus this method gives too large an estimate for the number 

of distinguishable configurations. With this assumption, the 

number of configurations available to the (n-f l)th polymer 
molecule is ^ v . 

= iN—mYz{z—lY-ya^N^-\ (3.51) 

where is the symmetry number of the r-mer molecules. It 
follows that the total number of distinguishable configurations 
of the molecules of the mixture is 


N/ n 

^ n=i 


(3.52) 


5129 


B the lactonal term is introduced to eliminate configura- 
which differ only in the interchange of two or more polder 
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molecules. From (3.51) and (3.52) we get for the niimber of 
configurations 




f2(2-l)’'-V) 

-Vr 1 1 

f (^V,>)! 1 

1 J 

.Y,! 1 

i(.Y/r-.^;):j 


(3.53) 


so that, to this degree of approximation, the contribution of the 
configurational terms to the free energy is 


ifT.Y.log- 




Y -L r Y 


feT.Y.log- 


N. 




+ 


-^kTX,log{z{z-iy-^!<r,}^kT{r-l)y\. (3.54) 

For the pure species of r-mer molecules there is, according to 
Flory’s method, a configurational contribution to the free 
energy amounting to 

feT.Y.log { 2 ( 2 - l)^-2;a,}-fer(r- l).A;+feTxY,log r. (3.55) 


firee 


kT 






y\+rN, 




rK 




(3.56) 


Introducing the volume fractions and v, of monomer and 
polymer respectively, this reduces to a particularly simple form: 


kT 


iN\logi-i—Y.logr,. 


(3.57, 


3.18. Huggins similarly computes the number of ways o; 
arranging successive pohnner molecules on a quasi-crystallin( 
array of sites. His analysis differs from that of Flory in that h( 
calculates approximately the effect, on the number of site: 
available for the segments of each polymer molecule, of tb 
blocking of sites by other, distant, segments of the same polyme: 
molecule. By considering the average occupation of neighbour 
ing sites by different pohuner molecules and by different seg 
ments of the same polymer molecule, Huggins finds that th 
number of sites available to accommodate the ?nth (2 < m < r 
segment of the 7ith polymer molecule is 

(2-l)(l-.^J 

l-2(r-l)i„/r2’ 

S.. = in—(77—1)7'(1 
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with <f>o a small constant which Huggins finds to be of the order 
of magnitude of 0*05. He thus obtains for the number of con¬ 
figurations of the nth polymer molecule 


-=Ui 


1=1 


-2(r-l)<f>,lrz 


With this value, he obtains for the free energy of nuxing, 


(Mk 

kT 


+ 




z 


(3.58) 


where z', which does not differ greatly from z, has been written for 


rz 




2 ^ 
1-^0 ■ 


(3.59) 


3.19. In the application of these formulae for the free energy 
of mixing to solutions of high polymers it must be remembered 
that the quantity r is very large, of the order of magnitude of 
many thousands. With this in mind, it can be seen that equations 
(3.47) and (3.58) are virtually identical and that equation (3.56) 
is the limit to which both tend as z approaches oo. In view of 
what has been noted in § 3.15, it would appear that the detailed 
consideration of the disposition on the array of sites of each 
segment of each polymer molecule is unnecessary and merely 
introduces unwanted algebraical complications.f Another 
determination of the combinatory factor by calculating the 
number of permutations of molecules on an array of sites has 
been given by Miinster.J He introduces the idea of a virtual 
molecide which consists of a group of monomer molecules of the 
same size and shape as each polymer molecule and considers the 


t A comparison of these three methods of determinim? 
fectop has hem given by Gee (Advances in CoOoid Science ,: 
a ^eral review has been given by Manster (KoU. Zeits 18 



5J KVAIAWTION OF THK ( OMHINATOKY FAiT'OH <>. 

iiunil>t'r ot juTmutation;^ ol jw>lyiiu'r inoUH.^uk*^ aii«l \ irtUAi molr- 
cules, allo\% iiijr for the num>K*r of wnfigurationi* of each Tnole<nile 
w hoo one ot its M*^inents is tixe<l. Ft>r th«* fre<" eiierffv of ililutHm 
his result is identical v^ith a first ap{»n»ximatioii oKtairn^l bv 
ex|mnding e<|uati<»n (3.47) as a ]>«»\%er s<*ri»-s in itK' xcdunte 
fraction of the {xjlyinc'r. 

A (juite different in<Klel. applicable t<.> a liquid near its critical 
fH>int. can also Ik* Usc*d ti> evaluate the <*ntropy of mixing.+ It 
is f'uuiu^l that, i»n this mmlel, th<* iinj_K»rtant i»aranieter is tl*e 
ratio of the molecular volumes, and if e<^uHtiim (3 47) is reuntt<-n 
in terms of the volume fractions of the coin|Kim*ntK. it is found 
to iilentical with the equation \%hich is obtaiiuxl from this 
mcxlel of a higlily im{>erfec‘t gas. 

It is noteworthy that determinations of the free energy of 
mixing on the V)asis of two quite ditTei\*nt mfKiels, one of uhich 
is a good representation of a liquid near its nieIting-|M>int and 
the other of which is a good representation of a liquid near its 
critical ]K)int, both lead to the same result. This increases the 
assurance with whicli equation (3.47) can l>e accepteil. Further¬ 
more, this result can now be deduced in a more rigorous way. 
To do this it is necessary to make use* of scune general statistical 
results which have l>een derived^ for simple chain, branched 
chain, and closed ring molecules. 

Basic statistical results 

3.20. The physical model which is considered is the quasi- 
crvstalline arrav of sites described in § 1.0, and the parameters 
wMch specify the composition of the mixture and the molecules 
of its constituent species have been defined in 2.2 to 2..>. For 
the present, attention is confined to mixtures of any numl>er of 
species for which the energy of mixing of each pair is zero, so 
that there is a random distribution of molecules on the array 
of sites. The total number of closest neighbour interactions in 
which molecules of species i are partners depends both on the 
number of molecules of this species and on the number of closest 

+ Hildebrand. J. H., J. Chem. Phys, 15 (1947), 225. 

+ Miller, A. R., Proc. Cambridyt Phil. Soc. 43 (1947), 422. 
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neighbour interactions which each molecule of this species yields. 
Thus the relative frequency of i—j interactions, for all values 
of » and j, is qiqfNiNj. The relative frequencies of various 
laical modes of occupation of a given pair of closest neighbour 
sites (conveniently referred to as id and B) are given in Table 2. 

I 

Table 2 

Rdative Frequencies of Typical Modes of Occupation of a 

Given Pair of Sites 



The definitions of and Xfj imply that the frequency of occupa¬ 
tion of the given pair of sites in the mode a is proportional to 
Xf^ and that the sum of the frequencies of occupation in the 
modes )3 and y is proportional to X^^. X^^ is proportional to the 
sum of the frequencies of occupation in the modes and y 
because these are the two possible ways in which the given pair 
of sites can yield an i—j pair. It therefore follows fit}m the 
equations (2.3) and Table 2 that 


X.- 




2 2 

k 

qig,N,N, 

k 


(all i). 


(3.60) 


(all i and j with » j) 


It will be noticed that the last line of Table 2, referring to the 

mode of occupation 8, follows immediately from equations (2.5) 
and (3,60). 

The probability that a given site, say A, is occupied by a mole¬ 
cule of species i is r^N^/2 AH the ways in which site A 

k 

can be occupied by a molecule of species t is included in the modes 
of occupation indicated by a, and 8 in Table 2. In this con- 
nexioi^ the mode of occupation j8 is regarded as typical of all 
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modes in which site A is occapied by a molecule of species i and 
site B is occupied by a molecule of any species other than i. 
Thus when we are considering all the ways in which site A r»a.n 
be occupied by a molecule of species t we have to consider the 
mode jS for all values of j not equal to t. It follows therefore that 
the probabilities of the occupation of a given pair of sites in the 
various modes indicated are K times the relative frequencies 
given in Table 2 where the constant K is given by the equation 

^ 2 9k^k 2 ^k^k — 

k k 

Without further elaboration it is therefore clear that the proba¬ 
bilities of occupation of a pair of sites in various typical ways for 
a mixture containing any number of species co nsisting of simple, 
simple chain, branched chain, or closed ring molecules are given 
in Table 3. 


Table 3 


Probabilities of Occupation of a Given Pair of Sites in 

Various Ways 



Mode of occupation 


Both by % (excluding adjacent groups of 
the same molecule ») 

A by i, B by/j 

A hyj, B by i/ 

By adjacent groups of the same molecule t 


Probability 


k t 

9i 9j NjlX 9k ^k 1 
(ri-gi)Nil2^k^k 

k 


It will be recalled that in these formulae the parameters are 
defined by the fact that zq^ is the number of closest neighbour 
interactions of a molecule of the species t, and that its value has 
to be determined according to the type of molecule. It should 
be emphasized again that these results follow from the physical 
meanings of the various parameters which have been introduced 
and do not depend upon the particular values which any of these 
parameters assume for any particular kind of molecule. We 
have already seen that for simple chain or branched chain mole¬ 
cules the parameters q^ have the values given by equation (2.1), 
while for closed ring molecules they have the values gj given by 

equation (2.2). 
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Statistical derivation of general formulae 

3.21. The basic differential equation can be derived by a 
statistical argument in the following way.f For simplicity, the 
argument will be given for a binary mixture. Consider an 
assembly which consists of a mixture of monomer and r-mer 
molecules in contact with their vapour. Define two states of this 
assembly in the following way: 

State a: monomer and r-mer molecules, 

state j8: N^-j-r monomer and N^—1 r-mer molecules. 

For a random mixture the relative frequencies of occurrence of 
states a and /3, denoted by and Pp respectively, are propor¬ 
tional simply to the number of distinct configurations of the 
ffiven numbers of molecules of each type. Thus 


P. 


g{Nu K) 


(3.61) 

Pfi g(I^l+r,JSr-l) 

an equation which merely makes use of the statistical weight. 
In the present case the randomness of the mixture imposes an 
additional constraint and the sum over all the states of the 
assembly reduces to the combinatory factor which is the degen¬ 
eracy of the state which has the given (constant) value of the 
configurational energy. 

If the total number of sites is kept fixed and is decreased 
by unity so that is necessarily increased by r, then the incre¬ 
ment in log g{Ni, N^) is given by 


[Alog9(i\^i,i^,)]^=_i 


l)-log g{Ni, K) 


log 


giNr, K) 


^(.S^i-l-r.^-l) 


TiaiTig 


8\ogg{Nj^,N^) 

dN, 


1 +log 

Jw 


p 


p 


0 . 


(3.62) 




Now consider the occupation of all but a selected group (F) 
of sites, r in number, which axe iuterrelated in such a way that 
they can accommodate an r-mer molecule. Let i 


onomer 


t Miller, A. R., unpublished results. 
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molecules and —1 r-mer molecules be arranged in any con¬ 
figuration, specified by C^, on the remainiiig N—r sites: denote 
the probability of the particular configuration by P,. The 
occupation of T by an r-mer molecule is directly proportional 
to Nf. Then, since there is complete randomness in the occupa¬ 
tion of the sites, the relative fiequency of the configuration in 
which an r-mer molecule occupies T and the remaining N—r 
sites are occupied in the manner C, is 


PrKIN, 


(3.63) 


apart fi:om a factor which may depend on r, q, z, and T but is 
independent of and N,. This is the firequency of occurrence 
of a state in which there are monomer and r-mer molecules 
in the complete assembly, that is, the state a. The probability 
of the occupation of F by r monomer molecules is, by an exten¬ 
sion of the result given in lines and y of Table 3, 


N 


\N^+qNj ’ 


(3.64) 


Thus the probability of the configuration in which monomer 
molecules occupy T and the remaining sites are occupied in the 
manner C^is AT / xr \r-i 

This is the probability of a state in which there are J^i+r mono¬ 
mer molecules and N^— 1 r-mer molecules, that is, the state /S. 
Thus the ratio of P„ to Pp can be obtained from equations (3.63) 
and (3.65); the result is independent of i^, that is, it is inde¬ 
pendent of the configuration of molecules on the sites outside F, 
and so must be valid whatever that configuration. Thus 


P. 


P 


P 


K 
iS' 


uN^y-^ 


(3.66) 


where p may depend on r, q, z, and T but is independent of 
and Nf. Substitution of the value of PJPp given by equation 
(3.66) in equation (3.62) yields the differential equation 


8 log g(Ni,N^} 

8 N. 



log p-f r log(i\r-ri\^,)—logi^^. 


(r — 1 )log{J\r— (r-? )i^,}, (3.67) 
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where has been replaced by its value in terms of 3* and 2^, 
which are taken as the variables. In integrating equation (3.67) 
the constant of integration is determined by the condition that 
when is zero, g{N^,N^) must reduce to unity. Integration of 
equation (3.67) yields, for the combinatorj’ factor, the expression 

logg(N^,N,} = _logxYi!-logi\;! +—log(.Vi+3-V,)-'- 

r—q 

_ ^ log( a;+ r. a; )! 4- log . (3.68) 

r—q 

In this formula the parameters q have to be given the values q' 
or 5 ', defined by equations ( 2 . 1 ) and ( 2 . 2 ), which are appropriate 
to simple and branched chain and to closed ring molecules 
respectively. For simple chain and branched chain molecules 

z{r-q) = 2 (r-l), 


SO that 



(3.69) 


If the values of the coefficients of the terms in equation (3.68) 
are replaced by these values, then it is easily seen that equation 
(3.68) is identical wth equation (3.49) wth put equal to zero. 
Thus the derivation of equation (3.68) by the statistical method 
of this section leads to the same result as the mathematical 
extrapolation by which equation (3.54) was first constructed, j 

3.22* The physical significance of p can be determined^ bv 
considering the case in which > r.^;. In this case (r-q)X, 

can be neglected compared with X and equation (3.68) becomes 
asymptotically equivalent to 



Provided A\ > rA; the factor which multiplies in the right- 
hand member of equation (3.70) is the number of ways of 
choosing a; out of N sitee. It may therefore be interpreted as 
the number of ways of choosing sites for one selected element in 
each of the N, r-mer molecules. Consequently the parameter p 
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must denote the number of distinct configurations of a poljuncr 
molecule -w-hen one of its elements is fixed on a given site. This 
verifies that p is a configurational factor which depends only on 
the geometry of the molecules and of the array of sites. It is also 
related to the factor/^ which was introduced in § 3.14 to allow 
for the contribution per polymer molecule to the combinatory 
factor due to the flexibility of the polymer molecules. The rela¬ 
tion between p and is given by 

0’rP = z/r- (3.71) 


Although it is not neeessary to evaluate p or in order to de¬ 
termine the thermodynamical properties of the mixture,f it is 
of interest to note the values assumed by p for a number of simple 
cases in which its value can be deduced from a knowledge of its 
physical significance. In the trivial case of a monomer molecule 
p is unity. For a dimer molecule p is zja, where a is the symmetry 
number of the molecule. Substituting this value of p and r — 2 
in equation (3.68) one recovers equation^ (3.28). In the case of 
trimer molecules three cases have to be considered. If the mole¬ 
cules are linear and rigid p is z'a: if thej^ are linear and flexible p 
is z( 2 —1)/<T. Substituting these values and r = 3 in equation 
(3.68) one recovers equations§ (3.41) and (3.41'). If the trimer 
molecules are non-linear and rigid, then p is zz ja, where z is the 
number of alternative sites for the third segment when the first 
two segments of the trimer molecule have been placed. For 
molecules which occupy more than three sites there are several 
alternatives. The minimum value of p, which occurs for linear 
rigid molecules, is z/u. An upper bound to possible values of p is 

2(2 _This value must be slightly greater than the precise 

value of p for flexible molecules because it includes certain con¬ 
figurations in which the chain molecule bends back on itself and 
two segments occupy the same site; such configurations should be 
excluded, but the error due to their inclusion is probably small. 


t Miller, A. R., Proc. Cambridge Phil. Soc. 39 (1943), 62-3; Guggenheim. 

E. A., Proc. Roy. Soc. A 183 (1944), 208-10; cf. § 3-15 above 

t Chang, T. S., Proc. Cambridge Phil. Soc. 35 (1939), 26o. equation (33), 

Miller, A. R., ibid. 38 (1942), 109, equation (4). 

8 Miller. A. R., ibid. 38 (1942), 109, equations (33) and (39). 
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Kinetic derivation of combinatory formulae 

3.23. The general formula given in equation (3.68) can also be 
derived by kinetic arguments using the principle of detailed 
balancing.f Correct equilibrium conditions can be determined 
by considering an arbitrary elementary process and its inverse, 
since in equilibrium any single process balances its inverse 
process. For the present application of the principle of detailed 
balancing the elementary processes chosen for consideration are 

(i) the evaporation of an r-mer molecule fix)m a selected group of 
r sites and the simultaneous condensation on these sites of r 
monomer molecules from the gas phase, and its inverse process 

(ii) which consists of the evaporation of r monomer molecules 
from the same group of r sites and the simultaneous condensation 
on them of an r-mer molecule from the gas phase. The rate at 
which the first process takes place depends directly on the rth 
power of pi, the partial pressure of the monomer molecules in 
the gas phase, and the frequency of occupation of the group of 
sites r by an r-mer molecule. It is therefore directly proportional 


Pi 


N,-]-rN/ 


(3.72) 


The rate at which the second process takes place depends directly 
on the partial pressure p, of the r-mer molecules in the gas phase 
and on the frequency of occupation of F by r monomer molecules. 
Using equation (3.64) it is therefore directly proportional to 


Pr 


N, 


N^+rN, 




(3.73) 


In equilibrium these two rates must be equal, whence 


Pi 


K (N^+qNrY-^ 


(3.74) 


where C depends on r, q, and z and on the temperature T but is 
independent of and N^. 

The relation between the partial pressure p of a perfect gas 
and its absolute activity A is 


p = ^kTfiT), 


GuggeDhe] 
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where f{T) is the partition function of a molecule in the gas 
phase apart from the volume factor. Since A/p is a function of T 
only, equation (3.74) can be replaced by 

N^\ N, I ’ 

where p, which can depend only on r, q, and z, is independent 
of T and of and N/, and and respectively are used to 
denote the partition functions of monomer and r-mer molecules 
attached to specified sites. 

The Gibbs function is related to the absolute activities by the 
equation 

dG = p^dNi+pydN^ = fe!riogAi(i^^i+fe!riogA,.d^. (3,76) 




Thus, taking N and as the variables, 



(3.77) 


The combinatory factor is related to the thermodynamical 
functions by the equation 


A 

kT 



log g(Ny, N^)+Ni log <f>i-\rK log 4>r> 


(3.78) 


a-«8s nming the pressure to be small. Equations (3.75), (3.77), and 

(3.78) yield the differential equation (3.67) and the result already 
given in equation (3.68). 

3.24. Since the methods which have been used in§§ 3.21 and 
3.23 are concerned with systems with zero energies of mixing, 
the assumption of complete randomness of mixin g is vahd and 
does not introduce any error. There is, however, a concealed 
error in the interpretation of complete randomness in the neigh¬ 
bourhood of a site occupied in a specified way, and this implies 
that equation (3.64) is not exact. This concealed error is due to 
the fact that segments from different chain molecules may be 
competing for the same neighbouring site. It is, however, likely 
to be small. The assumption of complete randomness, though 
probably a good approximation for simple chain molecules, is 
not quite so good for branched chain molecules, and is least 
satisfactory for closed ring molecules. It will appear below that 
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appHed to simple chain molecules it involves the same sort of 
approximation as arises in the extrapolation (§ 3.14) by which 
the results for mix tures containing dimer and trimer molecides 
have been generalized to cover the case of mixtures containing 

polymer molecules. 

3.25. The methods which have been described in §§ 3.21 and 
3.23 can be extended immediately to examine systems which 
consist of any number of species containing any types of mole¬ 
cules. In fact, from the form of equation (3.68) and the inde¬ 
pendence, for simple chain and branched chain molecules, of the 
coefficients of the terms in this equation from r and g, the 
extension is obvious.f We obtain 

logff(A'i) = i2log(2?i-®^<)-— 

+ ( 3 - 79 ) 

i 

For closed ring molecules, on the other hand, 

2 (r—g) = 2r, 


so that 


r—1 ir—1 
r—g~2 r ^ 


r —g 2 r 


and in this case the coefficients of the terms in equation (3.68) 
are not independent of r and g. However, in all applications to 
high polymers, r ^ 1 so that (r— l)/r differs only insignificantly 
from unity; neglecting the difference between (r— l)/r and unity 
cannot cause as great an error as is inherent in the interpretation 
of complete randomness. Thus to the approximation inherent 
in the derivation of equation (3.68), equation (3.79) is valid also 
if some of the species consist of closed ring molecules. This 
general formula can be derived by a simple extension of the 
argument of either §§ 3.21 or 3.23, and such a mathematical 
treatment confirms that equation (3.79) is directly applicable 
to systems of any number of species which consist of simple 
chain or branched chain molecules, and, provided r is so great 
that (r—l)/r differs negligibly from unity, it is applicable also 
when some of the species consist of closed ring molecules. As 


t Guggenheim, E. A., Proe. Boy, Soc. A 183 (1944), 303, § 8. 
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we shall be concerned to apply these formulae only to binary 
monomer-polymer systems, the mathematical details for the 
general caset of many components will not be given here. It 
should be noted, however, that with our interpretation of com¬ 
plete randomnesst equations (3.54) and (3.68) for a monomer- 

polymer system apply in all cases with the appropriate values 
of the parameters r and g. 

t Miller, A. R„ unpublished results. J See § 3.24 above. 
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COMPARISON OF THEORY WITH EXPERIMENT 

4.1. We have seen (§ 1.14) that the partial potentials in the 
mixture (liquid phase) are related to those of the corresponding 
pure liquids by equations of the form 


/4 


, (8^G 



(all i), 


(1.24) 


where AG is the Gibbs free energy of mixing, which is defined by 


the equation 


G 


1 Ot+AG, 

i 


where the Gf are the free energies of the respective species and 
G is that of the mixture. The partial molar free energy of mixing 
■of the mixture by the species i is defined by 


AG.. = L 


’8A 


8K 



(4.1) 


where L is Avogadro’s number. 

We have also seen that the partial potentials in the vapour 
phase for the solution and the pure species are related by 
equations 

logli. (alii) (1.28) 

Pi.O 

provided the vapour behaves as a perfect gas. For equilibrium 
it is necessary that the partial potentials in the liquid and 
vapour phases should be equal. Thus the preceding equations 

(all*)- (4.2) 


Pi.o 


It loUows that the free energy of dilution can be calculated from 
experimental measurements of vapour pressure or of any other 
property which is related to it thermodynamically. Such quan¬ 
tities are the osmotic and swelling pressures and the depression of 
the freezing-point. Measurements of all of these have been made 
in studying solutions of high polymers and gels. Vapour and 
osmotic pressure measurements have nroved moat 
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4.2. Many determinations have been made of the vapour 
pressure of solutions of high polymers. The quantities to be 
measured experimentally are the vapour pressure and the con¬ 
centration of the solution. The latter can be measured either by 
analysis of the solution or by transferring a measured amount 
of liquid from a reservoir on to the polymer. The former was the 
method adopted by Stamberger,t who first pointed out clearly 
the precautions which are essential in making vapour-pressure 
determinations on such systems. The second method is more 
convenient experimentally and makes it possible to use a simpler 
type of apparatus. J 

Osmotic-pressure measurements on solutions can be extended 
to very high dilution. In this way measurements are obtamed 
in a region of the concentration range not otherwise accessible. 
In making measurements of the osmotic pressure§ both static 
and dvnamic methods have been used. In the former methods, 
time is allowed for osmotic equihbrium to be attained and the 
osmotic pressure is measured by the difference in hydrostatic 
pressure between cells containing the solution and the pure 
solvent. In the latter methods, the rate of osmosis is measured 
as a function of an applied hydrostatic pressure and the osmotic 
pressure is determined as the hydrostatic pressure at which 

osmosis just ceases. 

Swelling pressure does not differ thermodynamically from 
osmotic pressure but may be distinguished by the fact that 
the high polvmer is not dispersed, so that the experimental 
techniques are different. It provides a possible means of 
carrying out experiments in the region of small vapour- 
pressure lowering and has been used to obtain measurements 

on rubber. II 

The depression of the freezing-point of low molecular polymers 


t Stamberger, P., J. Soc. (1929), 2318. « 147 

+ Gee G and Treloar, L. R. G.. Traw. Faraday Soc. 38 (1942) 147. 

8 Slu G. V.. Z. physikal. Chem. A 176 (1936), 317; Ca^r, S. R., and 

Biortt R., / (IMS), 660, G... G T,,^. 

(1940). 1162; Mever, K. H., and Boissonnas, C. G.. Hdv Chim A^, 23 (19^. 
430; Flor^^ P. J.'. J. Amer. Chem. Soc. 65 (1943), 375; Fuoss, R. M., and Mead, 

D. J., J. Pfiys. Che77i. 47 (1943). 59. 

i Posniak, G., Koll.-BeOufU, 3 (1912), 417. 
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has also been measured,f but for dilute solutions of rubber in 
benzene the freezing-point does not differ perceptibly from that 
of pure benzene. Measurements of the freezing-point of rubber 
in camphor have been attempted,! but it has been 8hown§ that 
unless the experiments are carried out in vactio the rubber under¬ 
goes severe 

made in vacuo on solutions sufficiently dilute to be fluid the 
freezing-point of the solution is indistinguishable from that of 
pure camphor. 

From this brief outline it appears that vapour-pressure 
measurements and osmotic-pressure determinations have pro¬ 
vided the most useful data in the study of the thermodynamical 
properties of solutions of high polymers. Details of the apparatus 
and of the methods used can be found in a review by Gree.|| We 
now consider some typical results and compare them with the 
theory which has been developed in Chapter III. 


degradation and that when the measurements are 


Vapour-pressure equations 

4.3. For a mixture with zero energy of mixing and containing 
Nf monomer and r-mer molecules the partition function is 

g{N„ 

where g(Ni, N^) is the combinatory factor determined by equation 

(3.64) with set equal to zero and are the partition 

functions, apart from the configurational factors, of monomer 

and r-mer molecules respectively. Assuming the pressure to be 

small, the Helmholtz and Gibbs free energies of the mixture are 
given by 


A 

kT 



log g(N^, log ^i-t- V, log <f>,. 



T Itemp, A. R., and Peters. R. H., Itid. Eng. Chem 
(1942), 1192. 


Pu^rer, R. et al, Ber. 60 (1927), 2167; ibid. 62 (1929), 2628; Kautschui 


5 (1929), 129. 


§ Gee, G., Trans. Faraday Soe. 38 (1942), 109. 

II Gee, G., Adtianees in Colloid Science, 2 (1946), 145 

S12S 
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« 

Ihe &ee energy of mi^Tig is therefore given by 


AG 

kT 


= log g{0, N ;)+log g(Ni, 0) -log g{Ni, N^) 

It therefore follows, by differentiating equation (4.4) with 
respect to and N^, that the molar free energy of dilution of 
the solution by the solvent is 

= logN^-\zlog{N^-\-qNAMlz-i)log(N^-^rN,), (4.5) 
R.1. 

and that the molar free energy of concentration of the solution 
by the solute is 


AG, 

RT 


(i*- l)f log^i?^' - Izq log^i+«^' 


Nr 


Nr 


(4.6) 


It has been seen, equation (4.2) above, that, for each species 
the vapour-pressure equations are given by 


AGi 


RT log^ 

Tifi 


(4.2) 


Using the values of the partial molar free energies of mixing 
given by equations (4.5) and (4.6) in equation (4.2), we obtain 

the vapoiir-pressure equations 




(4.7) 

jPw 

" (Ni+qN,)*^ ’ 


Ir.= 


(4.8) 

Prfl 

■ {N^+qN,)*^ 


in terms of the mole fractions of the components.t lie depar 
ture of the vapour pressure given by these equations fion 
Raoult’s law is iUustrated by the curves plotted in Fig. 2, u 
which the relative vapour pressures of both solute and solv^ 
are plotted as functions of the molar fraction of the solute 

♦ MQler A R . Proc. Cambridge PhU. Soe. 39 (1943). 54, (29) ^ 

( 30 );^Sel^ E. A.^Roy- Sac. A 183 (1944). 203. eqaatK«.3 (6.1 

/ft 9.V /fi.5). and (6.6). 
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various values of r, the number of segments in the polymer. In 
this diagram the full lines give the relative vapour pressure of 
the solvent and the broken lines that of the solute. For each com¬ 
ponent of the solution the straight lines give the vapour pressure 
for a mixture which obeys Raoult’s law. The curves showing 

10 





0 0-2 0-4 0-6 0-8 

Molar fraction of polymer 

Eio. 2. The effect of differences in the sizes of the constituent molecules on the 
vapour-pressure equations of solutions. The full lines give the relative vapour 
pressure for the solvent and the broken lines for the solute. In e^ C^ 
curves fa) are for a mixture of dimers and monomers, curves (6) for a mixture 

with moo Momra, c^es (e) for a mixture of polymers often segments 

onomers. The strz^t lines are for a mixture which obeys Raoul^law. 

departures from this law are, in order of increasing curvature 
for mixtures of monomers with dimers, trimers, and polymer^ 
containing ten segments. These curves show the influence of 
^e effects on the vapour-pressure equations, and they indicate 

that for a value of r as small as 10 there is already a considerable 

departure from the laws of ideal solutions even when the energy 
01 mixing is zero. 

Let us introdooe v, to specify the volume ftaetion of the polv- 

mer (,-mer meleoules), so that 1 - e, is the volume faction of the 
iUoiioiner. Then 

tN^ 


V. 
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and, since, for simple or branched chain molecules 






(Jz-l)r = izg-l, 

% 

equations (4.7) and (4.8) can be written in the form 

j>i.o 


(4.9) 


2l. 

Pr.O 


{l-(f-g)v,/r}**«’ 


(4.10) 


which is appropriate to simple or branched chain molecules. 

The osmotic pressure can be immediately deduced. It is 
related to the partial vapour pressure of the solvent in solution 
and the vapour pressure of the pure solvent by an equation which 

can be written , 

J Vidir = ll3’log{Pi.o/l>i}. 

0 

where is the partial molw volume of the solvent (monomer), 
TT is the osmotic pressure, and the other symbols have the 
meanings already specified in earlier sections of this chapter. 
It has been shownf that over the entire range of concentration 
the volume of mixing for mixtures of rubber and benzene is 
accurately zero. It foUows that the partial molar volume of 
benzene in such mixtures is equal to the molar volume of pure 
liquid benzene. We therefore assume that is independent of 
__and has the constant value »,o, the molar 


volume of the pure solvent. This 


«1,0 Pi 


(4.11) 


t Poanjak. G., KoU.-Beihefte, 3 (1912), «7, 

L which the volunxe fracti^ of^ 

almost identical, so that these measurements cover the whole mngv 
centrstions. 
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for the osmotic pressure. For dilute solutions this reduces to the 
weU-known formula 


•n 


RTl^p^ 
^1,0 Pija 


here Apj is the vapour-pressure lowering. 

From equations (4.9) and (4.11) the osmotic pressure 
lution of a high polymer is given by 


TT 


- RT, RT, 

iz——log 


V 


w 


Ni-\-rN, V 


1.0 


N^+rN/ 


(4.12) 


In the region of osmotic interest, rN^ N^. Expanding the 
logarithm in equation (4.12) and neglecting terms containing 
higher than the second power of rN^KNi+rN^), yields 


TT 




' the volume fraction of the polymer 
written 


TT 


l^^\2 z/wj 


(4.14) 


where V, = rvi „ and is the partial molar volume of the polymer 
in solution. K c gm.cm.-® is the concentration of the solution 
and if, are the density and the molecular weight of the poly¬ 
mer, this equation mav be written in the form 


IT 


RT 

Mr 


+ 


N) 


RT 


c. 


(4.15) 


polystyrene 


4.4. Measurements of the vapour pressure of solutions of 
polystyrene in various solvents have been carried outf and the 
results for non-polar solvents can be compared immediately 
with the vapour-pressure equation given abovej as equation 
(4.9). The experimental points which were obtained for solutions 
of polystyrene in toluene and in benzene are plotted in Pigs. 3 
and 4. Both sets of measurements were carried out at 20“ C., 

t Baaghan, E, C., Trans, Faraday 5oc., 44 (1948), 495* 
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hose with benzene as solvent using polystyrene of molecular 
weight 500,000 and those with toluene as solvent using poly- 
tyrene of molecular weight 200,000; the molecular weight was 
stimated from viscosity measurements at various concentra- 
ions. The full line drawn in each figure is calculated from 



Fig. 3. Variation of the vapour pressure of a polystyrene- 
toluene mixture. X experimental points, vapour pressure 
increasing; + experimental points, vapour preesore decreas¬ 
ing (after Baughan); -theoretical curve (for z = 6); 

--system obeying Raoult*s law. 


equation (4.9), and so shows the effect of different sizes and 
shapes of the molecules according to the theory developed in 
Chapter IH without allowing for the interactioiis between the 
molecules. The dotted straight line is the relation between the 
vapour pressure of the solvent and the volume ftaction for a 
svatem which obevs Raoult’s law. The significance of these 
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results will be discussed in § 4.6 below after similar results foi 
rubber-benzene systems hare been presented. 



Volume fraction of benzene 


Variation of the vapour pressure of a poly styren e* 
mixture; X experimental points (after Baughan); 
beoretical curve (for z = 6)j — — —■ system obeying 

Baoult’s law. 


Measurements on rubber 

4*5. Similar results have been found for vapour-pressure 
measurements on rubber-benzene systems. 

With volume fractions of benzene less than 0.015 a Macleod 
gauge was used to measure the relative vapour pressures. For 
volume fractions of benzene up to 0-8, a region which had been 
covered in earlier measurements,{the relative vapour pressures 

t Gee, G., and Treloar, L. R. G., Trans. Faraday Soc. 38 (1942) 147 

(1^29). 2318; W. J.. Bee. Trai Chlm. 51 
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»f the solution and pure benzene were measured by a direct 
nanometric method. For solutions with volume firactions of 
)enzene exceeding 0-8 the direct manometric method was 
napplicable and values were obtained by allowing a rubber 
iolution to come into equilibrium with a solution of triolein in 



Volume fraction of benzene 


Fig. 5. Variation of the vapour pressure of the rubber-bcaizaiio 

system; X experimental points;-theoretical curve (for z = 6); 

__ sy s te m obeying Raoult’s law. 

benzene of known concentration. In the equilibrium state the 
vapour pressures of the two solutions are equal, so that the 
relative vapour pressure of the rubber-benzene system could be 
calculated from that of the triolein. These measurements 
extended over the range 0-65 to 0-95 of volume fractions rf 
benzene. The rest of the concentration range was covered 
by osmotic-pressure measurements. The results of all 
measurements are shown plotted in Mg. 5. Again the fhU hue 
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shows the vapour-pressure curve calculated from equation (4.9) 
and the dotted curve is that corresponding to Raoult’s law. 

Discussion of vapour-pressure curves 

4.6. These diagrams exhibit two features of note. The first 
is the wide divergence of the experimental results from 
Raoult’s law. The second is that the curve calculated from 
equation (4.9) is in fair agreement with the experimental re.sult.s, 
so that the greater part of the divergence from the laws of ideal 
solutions is accounted for by the theory' which has been given 
in the preceding chapter. It therefore appears that the differ¬ 
ences in the sizes and shapes of the molecules of the constituent 
species of the mixture play a predominant part in producing 
divergence from the behaviour of ideal solutions. Other sources 
of departure from the laws of ideal solutions will be considered 
in the closing chapter of this monograph, but it is clear that, 
although there may be some regions of concentration in w'hich 
they assume a greater importance, the effects due to them must 
rather appear as refinements to the theorj' which has been 
developed so far. It should, however, be noted that the measure¬ 
ments on polystyrene do not extend to the region of extreme 
dilution and that although the measurements on the rubber- 
benzene S 5 ^stem have been extended over the whole concentra¬ 
tion range, a plot of the relative vapour pressures does not 
provide the most stringent test in the region of very dilute solu¬ 
tions. We now present some alternative ways of comparing 
the experimental results with the theory. 

The entropy of mixing 

4.7. The experimental results in which the relative pre.ssures 
are measured make it possible to determine the partial molar free 
energies of mixing directly. To calculate the entropy of dilution, 
or the partial molar entropy of mixing, it is necessary to deter¬ 
mine the heat of mixing. If Ml is the heat absorbed on mix ing 

molecules of solvent with molecules of polymer, then the 
heat of dilution or the partial molar heat of mixing is given by 



(4.16) 



• A. 
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The entropy of dilution is therefore given by 

TAiSi = AHi—AGj. (4,17) 

where AG^^ is determined by equation (4.5). -Now, in terms of 
the volume fraction of the r-mer molecules, equation (4.6) 
can be rewritten as 

^ = log(l-t;,)-Jzlogjl—(4.18) 

Measiirements of heats of mixing of very low molecular homo- 
logues of rubber with various solvents plotted s^ainst the volume 
fraction of one of the components give an almost parabolic 
curve,t which is the form expected on Hildebrand’s theory. Thus 

ah = acv^Vy, (4.19) 

where a is a constant whose value according to Hildebrand:^ is 

given by ^ _ (VBi—VB,)*, (4.20) 


where and E,. are the cohesive energy densities of the two 
species. Equation (4.19) leads on differentiation to the result 

= (4.21) 

where is the partial molar volume of the liquid in the solution, 

so tl»« ^ ^ ^ 

In a similar fashion one can obtain a relation between the en¬ 
tropy, fi:ee energy, and heat of concentration; thus 

(4.23) 

17? T 


(4.22) 


^Sr Ag, 

IT* 


(4.23) 


• a system with zero energy of mixing the free energy 
,nd the entropy of mixing are related by the equation 

AG = - TAS, 

the partial molar entropy of mixi n g the theory gives 
ASi _ i„„/, 14.24) 


^ = log(l—V,)—izlogjl 


(4.24) 


t Ferry, J., Gee. G., and Treloar, L. R. G,. Trxtns. Faraday Soc. 41 (1945), 

Hildebrand, J. H., and Wood. S. E., J. Chem, Phya, 1 (1933), S17; 
brand, J. H., Solubility of Non-Eleetrolytes, Chemical Catalogue Co„ New Yorfc, 

1936. 
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corresponding to equation (4.18) for the free energy. A test of the 
theory can therefore be provided by plotting the sum of AjS^-ir, 
determined from equation (4.24), and AG^Tv^, determined from 
the vapour-pressure measurements, against the volume fraction 
of the solute. According to equation (4.22) this quantity should 
be constant. It should be noted, however, that such a test 



Fig, 6. aV^/RT for different solvents as a function of the volume 
fraction of rubber; x benzene, -f- toluene, O carbon disulphide, 

O chloroform. 


depends not only on the validity of the theory of solutions of 
high polymers which has been developed so far, but also on that 
of Hildebrand’s theory of the heat of mixing of Kquids. The free 
energies of dilution of natural rubber by a number of liquids 
have been calculatedf from vapour-pressure and from osmotic- 
pressure measurements, J and the entropies of dilution have 
been calculated from equation (4.24) for z = 5. The results are 
plotted in Fig. 6, from which it is clear that, for mixtures of 

t Gee, G., Advances in Colloid Science, 2 (1946), 165. 

r i the vapour-pressure data are: benzene (Gee. G., and Treloar 

L R. G TransJaroday Soc. 38 (1942), 147); carbon tetrachloride, chloroform 

^ (Boissonnas, C. G., and van 

ier Wyk, A. J.. Helv. Chxm. Acta, 20 (1937), 779); carbon disulphide (Stamber- 

L ’ osmotic-pressure data are- 

r912)%n)™ (Posnjak. G., Koll.-Beihefte, 3 
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natural rubber in a number of solvents, is nearly constant 

for a wide range of values of An altemalave method of testing 
the constancy of 




Tf^ 


has been used extensively by Huggins.t He has shown that a 
wide range of data on polymer-liquid systems can be represented 
by an equation of the form 


RT 






(4.25) 


where is the limit to which given by equation (4.24), 

tends as z approaches oo, that is. 


(AS^) 


CO 




(4.26) 


and V is assumed to be a constant. Equations (4,25) and (4.26) 
imply that if a function O is defined by 


then 




w 


>5 


furthermore, firom equation (4. 

RT 

region 


(4.27) 


good approximatioa to A^i- It is evident by comparing eqoa,- 
tions (4.24) and (4.26) that v win neither satisfy equation (4.27) 
nor will it be independent of In fact, equation (4.26) provides 
only a rather insensitive test of the thermodynamical data on 
polymer-liquid systems. This insensitiveness of the method of 
plotting specified by equation (4.25) can be d^onstrated by 
plotting d) as a function of the volume fraction of polymer. 
From such curves it is evident^ that neither the value of z nor rf 
aVJT could be deduced even approximately in this way. A 

t HugginB. M. L.. ^nn. Ncu, Tort Aeod. Sci. 43 (1942). l .Ind.Bng. Okmm. 

G.. « CcOaid Sa^ 2 (1946), 145. Kg. 11- 
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more direct way of testing equation (4.22) is by using measured 
values of the free energy of dilution with measured values of the 
heat of dilution to determine the entropy of dilution in accor¬ 
dance with equation (4.17) and to compare the results with the 
values calculated from equation (4.24). This we now do. 

4.9. Measurements have been madef of the heats of dilution 
of natural rubber with a number of liquids. For most of them 
AHJvl was sensibly constant throughout the w^hole range of 
volume concentrations over which measurements were made. 
These measurements were combined with the free energy of 
dilution of natural rubber by the same liquids, calculated from 
the available vapour-pressure data. The results for the hydro¬ 
carbons, in the hquids which were investigated (benzene, 
heptane, and toluene), are plotted in Fig. 7. The curves are 
calculated from equation (4.24) for z = 4andz = 6, It is evident 
that the estimated entropies of dilution are generally of the form 
predicted by the theory but are not in complete quantitative 
agreement with it. At least, however, these results indicate a 
semi-quantitive agreement between the theory of polymer 
solutions and the thermodynamical quantities determined 
experimentally. 

It should also be noted that some part of the disagreement 
may arise from errors in the heats of dilution determined experi¬ 
mentally. In fact, the heat of dilution of rubber by benzene, 
which was used in calculating the experimental points in Fig. 7, 
did not agree with an earlier determination by the same 
workers.f Further measurements were therefore made§ of the 
heat of dilution of natural rubber by benzene.|| There is a con¬ 
siderable scatter of the experimental points, so that the measure¬ 
ments are not accxirate enough to determine the precise shape of 
the curve. Nevertheless, they suggest that the heat of dilution 

t Ferry, J., Gree, G., and Treloar, L. R. G., Trans. Faraday Soc. 41 (1945), 
340. 

t Gee, G., ibid. 38 (1942), 418. 

§ Gee, G., and Orr, W. J. C., ibid. 42 (1946), 507. 

II Of the solvents for which experimental points are plotted in Fig. 7, benzene 

is the only one for which a re-determination of the heat of dilution has been 
carried out. 
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of natural rubber by benzene increases continuously with the 
volume fraction of the rubber, that is, decreases with increasing 



Flo. 7. Entropy of dilution of natural rubber with difierentaolvenfa. 

The full curves are from the theory with the values of s 

imliijoted. The experimental points refer to beniene (X) at 25°C>, 

heptane (O) at 40*C., toluene (+) at 26*C. 

dilution. From a reasonable curve drawn through the experi¬ 
mental pointef Gee and Orr recalculated the experiments 

4 Sm Fiff. fi bfiliSW. 
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value of the partial molar entropy. Tliis is represented by the 
full curve in Fig. 8, while the dotted curve is calculated from 
equation (4.24) with z — 6. It can be seen that although for 
volume fractions of rubber greater than about 0*25 the two 
curves have the same general shape, there is not complete agi^ee- 



Fig. 8. Variation of the partial molar entropy of dilution of 
natural rubber with benzene with the volume fraction of rubber; 
smoothed experimental results,-theoretical curve. 

ment between them. At low rubber concentrations there is an 
appreciable discrepancy between the shapes of the experimental 
and the theoretical curves. In particular, the point of inflexion 
in the experimental eur\ e in the neighbourhood of ?; = 0-25 
has no counterpart in the theoretical curve. This is the main 
point of disagreement between tlie theory and the experimental 
results; apart from tliis both cuiwes are of the same general 
shape. The agreement between the measurements of vapour 
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jressure and the vapour-pressure curve calculated from equation 
4.9) is evidence that the free energy of dilution accords reason¬ 
ably well with that calculated from the theory. The discrepancy 
between the entropy of dilution calculated finm the measured 
5 -apour pressure and that given by the theory must therefore be 


403 - 
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Volume fraction of rubber 


1-0 


no 9 Experimental values of the heat of dilution of nat^ ^ 

ilG. y. icnni^tic data. A osmotac data, 

benzene; X vapour-pressure data, -f- isopiestic oaia, v 

□ low vapour-pressure data. 

due to the measured fall in the value of Aff j/if as t; decrea^ 
as shown in Fig. 9. If a direct plot of AGJTtf agam^ v is m e 
from the measured vapour pressures there is no pomt of mflexi 
in the neighbourhood of iv = 0-25 (compare Figs. 4 and o above). 
In view of this and of the evident uncertainty m the measure¬ 
ments of the heat of mining, it is apparent that the 

“of diluXn. a. sho^ in Fig. 8. ahonld not be acconied 
too great weight. 

. .__fViPre is a further point to be noted. 
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therefore gives a value for the free energ}.* directly. In the 
experiments, measurements are made of the vapour pressure of 
the solution compared with that of the pure solvent. This 
measurement determines the free energy of mixing directly. 
With measurements of the heat of mixing, the entropy of dilution 
as plotted in Fig. 8 can be derived from it. It must, however, 
be noted, as is clear from the preceding remarks, that the 
experimental and theoretical values of the entropy of dilution 
plotted in Fig. 8 are not strictly comparable. Both in the experi¬ 
mental work and the theory the primary quantity is the free 
energy of mixing. Consequently, irrespective of the reliance 
which can be placed on the measured values of the heat of mixing, 
greater weight should be given to the agreement between the 
free energy determined by experiment and calculated from 
the theory (Figs. 3, 4, and 5) than to the divergences between 
the derived values of the entropy (Fig. 8). It is concluded that 
in view of the elementary nature of the physical model which 
is employed in the theory, and of the assumption of random 
mixing (which assumes zero energy’' of mixing), the agreement 
between the theory and experiment is remarkably good. Cer¬ 
tainly, the semi-quantitative agreement which has been found 
indicates that the main cause of the wide departures from 
Baoult’s law shown by solutions of high pol^uners is to be sought 
in the effects of differences in molecular sizes and shapes. 

4.11. It is clear, however, that it is desirable to determine 
the effect of a non-zero energy of mixin g. For if the point of 
inflexion which is given by the experimental res\ilts, and arising 
from the experimental values of the heat of mixing, is a real 
effect, then it should also be manifest when the effect of a non¬ 
zero heat of mixing is examined theoretically. To do tliis it is 
necessary to evaluate the more general combinatoiy' factor 

This is the subject-matter of the next chapter. 


5129 
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THEORY OF MIXTURES WITH NON-ZERO ENERGY 

OF MIXING 

5.1. In an earlier chapter of this work the combinatory factor 
g(N^) which occurs in the partition function of a mixture with 
zero energy of mixin g and which contains simple chain, branched 
chain, or closed ring molecules has been evaluated. In order that 
the effect of a non-zero energy of mixing on the vapour-pressm'e 
equations and the other thermodynamical quantities can be 
determined it is necessary to evaluate the combinatory factor, 
denoted by g{Nf, X,-^), which occurs in the partition function for 
mixtures with non-zero energy of mixing (equation (1.34)). It 
is clear that Bethe’s method could be used, for it was precisely 
to take account of interactions between neighbouring particles 
that it was devised. In determining giN^) a set of partial differen¬ 
tial equations, one equation to every component of the mixture, 
corresponding to equation (3.38), is obtained, and this set of 
equations is sufficient to determine g(Nf) as a function of the N^, 
"When there is a non-zero energy of mixin g there are, in addition, 
a set of partial differential equations, as many equations as 
there are quantities corresponding to equation (3.39). The 
complete set of partial differential equations (3.38) and (3.39) 
are sufficient to determine giN^, ^ij) as a function of the and 
the X,-,-, provided that the conditions of integrability of the form 


It 


(all i). 

(5.1) 

^ 1 


(all h j), 

(5.2) 


are satisfied. In these equations e,-, ej are the parameters intro¬ 
duced in Bethe’s method corresponding to the occupation of a 
site by a molecule of the type specified by the suffix and has 
been written for exp(—where is the energy of 
mixing. This method was first usedf to determine the com- 

t Chane, T. S., Proc, Cambridge Phil. Soc, 35 (1939), 265. 
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binatory factor for mixtures of two species of monomer molecules 
and of monomer molecules with dimer molecules. However, 
the results can be derived in a more straightforward manner 
naing the methods of chemical kinetics; in this connexion this 
method has been called the condition of quasi-chemical equili- 
brium.f This equation was first put forward in the consideration 
of non-ideal solutions of a class, defined as regular solutions 
(§ 1.12). The method is, in fact, equivalent to Bethe’s approxi¬ 
mation, and we now proceed to give an account of it and to 
demonstrate their equivalence. 

Condition of quasi-chemical equilibrium 

5.2. The configurational factor in the partition function for a 
strictly regular solution has been derived in § 1.13. With suitable 
choice of the energy zero it is given, for a binary mixture, by 


n(T) 




/ ...J 


(6.3) 


It was also pointed out in § 1.15 that if a quantity is defined by 


■WgX^ilkT j j 


J... J 


(1.30) 

then the evaluation of the integral in equation (5.3) is reduced 
to the determination of 


^jworkers in the theory of solutions commonly assumed 
zXi^ is the average number of closest neighbour i~~j pairs, 
sufficiently good estimate of is equal to X riven bv 


X^ = iN,^X){N,-X), 


(5-4) 


where now, since we are dealing with binary mixtures, the 
suffixes to Z c£m be omitted. From equation (5.4), X is neces¬ 
sarily independent of T. Since zX is the average number of 
^ J paim, the numbers of i—i and j—j pairs respectively are 
^{N^—X) and Z), the factor one-half being introduced 

to take account of the fact that each %—i pair involves two mole¬ 
cules of type i. The physical meaning of equation (6.4) is there¬ 
fore that the numbers of i—i, i—J, and J—j pairs are determined 

t GuggezLheiin, E. A., Proc. Bov. Soc. A oaa 
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simply by the numbers of molecules of different tj-pes. It thus 
corresponds to the assumption of completely random miTing 
and can be strictly vahd only for systems with zero energy of 
m i x i n g. It is therefore a somewhat crude approximation for 
other cases. Since the integral on the left-hand side of equation 
(1.30) is the configurational term for a random mixture of the 
two species, denoted by giN^), the configurational free energy of 
the system must be _ 

kT log g{Ni)+Xu'. (5.5) 

If the combinatory factor for a mixture with non-zero energy of 
mixing be denoted, in our usual notation, by gr(iV,-,X), then the 
configurational free energy of mixing is 


kT\ogg(Nf,X)-{-Xw. 


(5.6) 


Since, in general, g{Nf,X) is necessarily different from g{Nf), it 
follows that X and X are necessarily different. It has been seen 
that to assume that X is given by equation (5.4) implies that 
there is completely random mixing so that the combmatory 
factor must reduce to p(iVf). Thus if zX given by equation (5.4) 
is assumed to be the average value of the number of i—j pairs, 
then 


X = X. 


(5.7) 


From equation (1.32) it follows that when X is independent of 
T then XjT can differ from XjT by, at most, an additive con¬ 
stant independent of T. The physical argument given above 
indicates that in this case the constant of integration which is 
introduced when equation (1.32) is integrated must be zero. It 
will be seen below (§ 5.6) that when y(X,-, X) is evaluated and the 
value of X given by equation (4) is substituted for X, then 
g{Ni,X) does in fact reduce to giN^), which confirms equation 

(5.7). 

5.3. It has been seen that equation (5.4) corresponds to the 
assumption of completely random mixing. An alternative 
formula which has been proposedf is 

(Ni-X){Nj-X) = (5.8) 

where X is now a function of T and there is no longer a simple 

t Guggenheim, E. A., Proc, Roy. Soc. A 148 (1935), 304. 
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relation between X and X, as is evident from equations (5.5) and 
(5.6), Equation (5.8) is clearly based on the equations of chemical 
kinetics, and, even if it is not exact, it should at least be an 
improvement on equation (5.4), If the symbol i~j is used to 
denote a pair of sites occupied by segments of two different 
molecules, one of type i and the other of type then the re¬ 
arrangement of the molecules which occurs when an i —{pair and 


^j—j pair are destroyed and two new i—j pairs are created can 
be represented symbolically by 


(»■-*)+ ij-j) ^ (5.9) 

Equation (5.9) recalls the formula for a chemical reaction; if 
it did represent a chemical reaction between molecules, the 
equilibrium condition would be 



arrangement of molecules is not that of a chemical reaction, it 
is not unreasonable to expect that the equilibrium condition 
between t—t, i—j, andj—j pairs should take a similar form. It 
was on this basis that equation (5.8) was proposed in investigat¬ 
ing the properties of regular solutions. When a relation such 


as equation (5.8) or equation (5.10) is applied to the numbers of 

neighbouring sites which are occupied in various ways, it can be 

conveniently referred to as a condition of ‘quasi-chemical’ 

equilibrium. It should be emphasizedf perhaps that although 

equations (5.8) and (5.10) have the same form as the law of mass 

action in chemical reactions, it is the relation of a given molecule 

to the surrounding structure, not its relation to another mole- 
cule, which is being considered. 

In applying the method of chemical kinetics, by introducing 
the condition of quasi-chemical equilibrium, to derive the pro- 
prties o^regul^ solutions, Guggenheim did not distinguish 
between X and X. That their identification was invalid when 

IS given by equation (5.8) was pointed out by RushbrookeJ 


I W. L.. and Williams. E. J., ibid. A 

I Bushbrooke, G. S., ibid. A 166 (1938), 296. 


145 (1934), 669. 



6 


:MIXTinRES ^^^TH XOX-ZERO ENERGY Ch. V. § 5.3 

cf. § 5.2), who also showed that in this case the condition of 
[uasi'Chemical equilibrium could be derived from the approxi- 
nate formula for the grand partition function obtained by 
Lethe's method. It was later shown that, for particles each of 
rhich occupies only one site of the array, the condition of quasi- 
:heniical equilibrium and Bethe's equations are equivalent for 
ill regular assembliesf and also for assemblies with long-range 
)rder. j More recently, their equivalence has been shown also for 
issemblies of particles each of which occupies more than one site 
)n the arrav.§ The latter demonstration, which is the one 

k 0 

mmediately releTant to the purpose of this monograph, will 
low be considered. 


Equivalence of Bethe’s method and the condition of 

quasi-chemical equilibrium 

5.4. The equivalence of Bethe's method and the condition of 
quasi-chemical equihbrium for mixtures which contain polymer 
molecules can be shoisii in the following way. Let (i, i, S) denote 
the term in the grand partition function, constructed according 
to Bethe's method, in which the central site of the group of sites 
which is considered in detail and a selected first shell neighbour 
are occupied by segments of different molecules of ty-pe i and 
all the remaining sites are occupied in some specified manner 
denoted by S: define {j.J.S) in a similar manner. Let 
denote the term in which the central site and the selected first 
shell site are occupied by segments of molecules of types i and j 
respectively with aU the other sites occupied in the maimer 5; 
likeviise let {j. i. S) denote the converse manner of occupation of 
the central site and the selected first shell site. One of the 
essential approximations of Bethe's method is the assumption 
that the ratios (i, i, S) (!,J, S) and (J, i, S) (jJ, S) are indepen¬ 
dent of S. This follows from the way in which the parameters ot 
Bethe’s method are introduced and is clear from the considera- 


+ Guggenheim. E. A., Proc. Roy. Soc^ 17^11940) 189 

* Fo%vler. R. H., and Guggenheim. E. A., ibid. A 174 (1940), 189. 

F A ibid. A 183 (1944), 213. 
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tions of § 3.7. In fact, the parameters e and tj which occur in 
Bethe’s method can be defined by 


(i, i, S) e 


(5.11) 


(j, i, S) 

{3,3, S) 


(5.12) 


A further set of quantities can be defined in the following way. 
Denote the sum of all the terms in the grand partition function 
constructed according to Bethe’s method, corresponding to the 
occupation of the central site and the selected first shell site by 
segments of different molecules of type i by {i, i ); define {j,j) in 
a like manner. Denote by {i,j) the sum of all the terms in the 
grand partition function in which the central site is occupied by 
a segment of a molecule of type i, and the selected first shell site 
is occupied by a segment of a molecule of type j ; hkewise define 
{j, i) as the sum of all the terms corresponding to the converse 
manner of occupation of the central site and the selected first 
shell site. Since the left-hand members of equations (5.11) and 

(5.12) are independent of S, then summation over aU configura¬ 
tions S yields the equations 


(i, i) ^ e 
{i,3) V 



(5.13) 

(5.14) 


The central site and the selected first shell site are chosen arbi¬ 


trarily; consequentlydhey are necessarily physically equivalent. 

Their physical equivalence implies 


and consequently 


(b 3) = (3, i), 
{{i,3)+{3,i)}^ 

Mi,i){j,j) 


(5.15) 

(5.16) 


The frequency with which a particular mode of occupation of 
the sites of the array occurs is proportional to the corresponding 
term in the grand partition function. Thus (i, i), which is the 
sum of a number of terms in the grand partition function, piust 
be proportional to the frequency with which the selected pair 
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of sites is occupied by segments of two ditferent moleculo of 
type i; likewise for (/.j). (J. {). and (j.j). The central site and the 
selected first shell site are chosen in an arbitrary manner and are 
physically indistinguishable from any other pair of sites in the 
array. Thus the frecpiency with wliich they are occuj.ied in a 
given mode is proportional to the frequency with whicli closest 
neighbour pairs, of the kind specified by that mode of occupation, 
occur throushout the assemblv. Thus 


= J,.,.: : AX . 


( 0 . 1 1 ) 


Ecjuation (5.17) identifies equation (5.16) \\ith the condition of 
quasi-chemical equilibrium provided only that 


which is precisely the value assigned to 7 ; in Bethe's approxima¬ 
tion. This establishes the equivalence of the two methods. 


Free energy of a polymer-solvent system 

5.5. We now consider how the condition of quasi-chemical 
equilibrium, that is. the method of chemical kinetics, can be 
applied to binaiy mixtures in which one of the constituents is a 
high polymer.T The generalization of the condition of Cjuasi- 
chemical equilibrium, equation (5.S), to a binaiy* mixture of 
monomer and r-mer molecules is 


(A\-A-)(9A;-A') = (5.19) 

where, as previously defined, ir is the energy- of mixing and zq is 
the number of closest neighbour interactions in which each r-mer 
molecule takes part. The solution of equation (5.19) is 


A* 


•2X* 


(5.20) 


where X* is the limit to wliich X tends as ic approaches zero. 
Thus A^=^ is the solution of the equation 

{X^-X){qX-X) = 

which is appropriate to completely random mixtures. This gives 



t Orr, W. J, C., Trans. Faraday Soc. 40 (1944), 320. 
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Also, j3, which is defined by 






(5.21) 


is a function of T. From equation (1.32), X is then given by 


l/T 


X = 2X*T 



d(llT) 
l+)3 


(5.22) 


CO 


Changing the variable of integration from l/T to ^ yields 


% 


4 



2X*zkT 


P 


w 





(1+^)(^2_^2)> 


where k is defined 




lN,-qNr V 

U+gATj • 


Integration of equation (5.23) gives the resultf 


(5.23) 


1 = £!^/_i_Iog_ 

w U+K ^(1-k)(14-^)^1_^ S( 1 _^)( 14 .^)J’ 

(5.24) 

where the constant of integration has been determined from the 

feet that, when the energy of mixing is zero (that is, jS is unity), 

X and X are equal, and each reduces to X*. That the value of 

X given by equation (5.24) satisfies the required conditions can 

be verified by using equation (1.32) to obtain X from it with the 
help of the relation 

W j32-/f2 

dil/T) ~zk ^ 

and thus recovering equation (5.20). This general result was 
given by Orr; an approximate solution for a mixture of two 
different kinds of monomers had been given earlier by Rush- 
brooke.J By expanding the radical in equation (5.21) and 

t Ibid. 

j Rufihbrooke, G. S., Proc. Boy. Soc. A 156 (1938), 296. 



90 


MIXTUBES XOX ZERO ENERGY Ch. V, | 5.5 

retaining the first three terms, an expression for i equivalent 
to Rushbrooke's approximation is obtained; thus 





(A\-gA;)(,d-l) 


1-/3 


2(e 


zhT 


1) 


I 




2'r'rk T 


1 ) 


i 


To the same order of approximation, we have 


(5.25) 


A = A 


A 


zkT 


yj_—qy,\ 2/f 


(f 


2’C zkT 


1)-1 


J 


(5.25') 


in which the constant of integration has been determined bv the 

_ % 

fact that for zero enersrv of mixing A reduces to A*v, 


5.6. From equation (1.30) the contribution to the free energy 
of mixing due to the non-zero energy* of mixing is given by 

AJ, = A(?, = Xh% (5.26) 

while the heat or energx’ of mixing is given by 

AH = AF = Aic. (5.27) 


In addition to AI^ or AG, arising from the non-zero energy of 
niixin<7. there are contributions to the free energy of mixin g due 
to the combinatory factor in the partition function; these are 
given by equation (3.62). The total configurational free energy 

of mixing: is therefore 



where A^"lc — ^ogglAj. A,.) g(0,A,). 

From the remarks made in connexion with equations (5.5) and 
(5.6) above it will be clear that this use of the condition of quasi- 
chemical equfiibrium is, in fact, equivalent to the evaluation of 
the combinatoix- factor glA’., A.-^-) in terms of the combinatory 
factor g(A;) which has been determined previously. Presentation 
of the explicit formulae to which equations (5.27) and (5.28) with 
equations (5.20) and (5.21). or approximately with equations 
(5.25) and (5.25'). lead will be deferred until the combinatory 

t Misprints in Orr's equations (10) and (14) should be noted, as a 
take in sign in the expression for N immediately foUowmg Kushbrooke = 

equation (6). 
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factor g(Ni, X^j) has been determined directly in terms of g(N{) 
with the help of the condition of quasi-chemical eqiuhbrium.f 

5.7. The evaluation of g(Ni,Xij) by direct methods is pro¬ 
hibitively difficult in the general case. However, it is possible 
to set down certain general considerations which g{N^,X^j) 
must satisfy, and, with the help of the conditions of quasi¬ 
chemical equilibrium, these are sufficient to determine the form 
of g{Ni, Xij) unambiguously in terms of g{,Ni). We now state the 
conditions. 

(i) To ensure that the free energy of a system of a given tem¬ 

perature and composition shall be proportional to the size of 
the system it is necessary that {log must be 

homogeneous and of zero degree in the 

(ii) To ensure that the free energy of the assembly is a con¬ 
tinuous function of its configmational energy W in the neigh¬ 
bourhood oiW = 0 it is necessary that, when W is zero, g{N^, 
should reduce to g(N^, apart from terms negligible compared 
with Ni. 

(iii) It will be assumed that the conditions for maximizing 


g{Ni, Xij)eTL^{-W{Xif)lkT} 

with respect to X^j are equivalent to the conditions of quasi¬ 
chemical equilibrium between the alternative methods of occu¬ 
pation of a pair of neighbouring sites. The argument can be 
presented generally for a mixture with any number of compon¬ 
ents of molecules of any type—simple, simple chain, branched 
chain, or closed ring. In presenting the general argument the 
notation which has been specified in Chapter II is used. 

The condition of quasi-chemical equilibrium, equation (5.10), 
can be written in the form 


(all i,j with i ^ j) (5.29) 
by using equation (2.3). Quantities Xf^ are introduced with 


t Gu^enheim, E. A., Proe. Roy. Soc. A 183 (1944), 213. 
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meanings con'esponcliiig to that assigned to iJ* in ^ 5.5. That 
is, X* is the solution of the equation 


k k 

A: 

and likewise the quantities X* are given by 




whence 


(5.30) 



(all 0. 


(5.30') 


It is now possible to construct an expression for the number of 
distincruishable conficnirations of all the molecules in which the 

V- 

confifrurational energy- has a given value, that is. the J,-,- have 
given values. The formula for g(X^, X^j) determined by the three 
considerations set out above is 


log^U\. J,/) = logguV,.)- 


viocriihil'^ 

^ ^ r - T r ‘ 


{zxfx'y^’i 


ij 


• ) 


(5.31) 


where {/(V,) is given by equation (3.6S). It can be seen from 
equations (5.30) that the are homogeneous and of the first 
degree in the iN'-, from which it follows that {logjliN’,-. J,j)} J 

given bv equation (5.31) is homogeneous and of zero degree in 
the Xi- From the definition of the X* it is emdent that when all 
the i/.-.j are zero reduces to g{Xi). Finally, if 


log{?(35.V,-,-)exp( 


IT kT)] = hsgiXi^Xij)- zkT 

iO' 


is maximized with respect to X,-;, the condition which is obtained 




— U'.j zkT = 0 {i —j), 


which is precisely the condition of quasi-chemical equiUbrium, 
equation (5.29). Thus equation (5.31) for ^(iN,-,X,-j) satisfies the 
conditions required for the combinatory factor. It can be 
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simplified by substitutiiig in equation (5.31) the values of 
given by equation (5.30) and of giN^). For simple or branched 
chain molecules equation (5.31) assumes the form 


= 2 ^i^ogpi— 2 iogiv;.!+z 2 iog(?<iv;)! 


(|2—l)log(2] riN^\-\z ^ log^j^iV^ 



'X 


ik 


zJ'logXf^. (5.32) 


A 


&ee energy of the assembly is given by 

-fcT 2 4>i-kT log g{Ni, X^^)+z (5.33) 

» i.i 


where is the partition function for a molecule of type 
kinds of energy apart from the configurational energy, 
denotes the equihbrium value of determined by e 
(5.29). Equations (5.32) and (5.33) yield, for the free 
of the assembly (at negligible pressures), 


kT 


4 t i { 

+(¥- i)iog(2 '+^2 2 Jog(?<-^<- 2' -^ i *)+ 


+z 2' log + 2' 

imM 


(5.34) 


From this all the thermodynamical properties can be derived. 
In particular the partial potentials and the absolute activities 
A* can be obtained by differentiating equation (5.34) with respect 
to iV,. Since SX/aX<y is zero for all i,j, variations in X^j may be 
ignored in carrying out this differentiation. Thus, one derive. 


kT 


logAf 


log tf>i Pi-^-logNf—zgf log 


-f (|z-l)rilog r^N^^^^\og^^Ni~ ( 5 . 36 ) 


which leads to 




je (I 


^iPi 


iaiN^fu 


for simple chain and branched chain molecules. 


(5.36) 
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These equations for the partial potentials lead immediately 
to the vapour-pressure equations. Specifying the partial vapour 
pressure of the component % in the mixture by p,* and the vapour 
pressure of the pure liquid by then 

^ _ r,. xA; 1 

= r (5 3 -) 


k 




In these formulae the values of the are given in equation 
(2.1) and the are given by the solutions of the equations of 
quasi-chemical equilibrium -which are determined by the set of 
equations (5.29). These cannot generally be solved exphcitly. 
Further simplification is possible only for a binary misture. 

5.9. When there are only two types of molecules, specified by 
1 and j in the misture, equation (5.32) reduces to 


g{N,,N^,X) 




X)!p 

(5.38) 


and equation (5.34) reduces to 


G 

kT 


A 

kT 


Xi log Pi - log pj ^ log JN;- ! -f log X ^! - 

2 log( 9 i Xi) ! -2 log(g,- iV;-)! -f (iz-1 )log(ri X^+r^ Xj )! -}- 

2 log X!+iz log(3i .A’i-X): -f i2 Iog(g; x^-xy.A- y^x, 

(5.39) 


where the suffixes have been dropped from X;, as there can 
be no risk of confusion by so doing. Formulae equivalent to 
equations (5.38) and (5.39) had already been obtained in th^ 
special cases when Guggenheim deduced the general formulae 
for mixtures containing simple or branched chain molecules. 
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Chaaigt had used Bethe’s approximation in a way similar to that 
described in § 3.7 to examine the two simple cases 

Ti = Tj — = 1 

and 

’•<"=2, qi = 2(z-l)lz, Pi = zl2, r^. = q^ ^ p^ ^ 1. 

OrrJ had used the condition of quasi-chemical equilibrium as 
described in §§ 5.5 and 5.6 to obtain the free energy of mixing for 
a system containing monomer and simple chain molecules, that 
is, for the case 

U = qi = {rz-2r+2)lz, Tj ^ q^ ^ p^ = 1. 


5.10. For this latter case the absolute activities are given by 
equation (5.36) with the appropriate values of the and j,-. 
Using the notation of § 5.5 for such a binary mixture we obtain 


Ax 


N, f(N,+rN,){N,-X)\i- 
<f>^{N,+rN,)\ Nl 


(5.40) 


where is the solution of equation (5.19). Using the parameters 

K and p, which have been introduced in § 5.5, these expressions 
reduce to 


- I 2(N,+rN,)(^-^) 

1 ^i(i\^i+ri\;)l(i\rj+giVr)(l+;3)(l-/f)/ ’ 





•l>r Pri^i+rK) \ (Ni-\-qN^){l +P){1 +k)I 


(5.43) 


^Aj^o, Q denote the values of the absolute activities for the pure 
substance we obtain 


A. = _^i / ^N^+rK){p-K) 

Ax.o N^-\-rN^ l(-A^i+3-^r)( 1 +)8)( 1 -ic)/ ’ 

I 2(N^+^i)W+k) m 

Kfi NJr-\-N, \{JffJq-{-N,)(l +;3)( 1+)/ ' 

t Chang, T, S., Proc, Cambridge PhU, Sac, 35 (1939), 265, 
t On*, W- J. C., Tram, Faraday Soc. 40 (1944), 320. 


(5.44) 

(5.45) 
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The relation between the partial potentials and the absolute 
activities is . 

^ = exp(fi,—o) kT, 


and. for equilibrium with the vapour phase, the right-hand side 
of this equation is the ratio of the partial pressure of species i 
above the mixture to the vapour pressure of the pure species i. 
Thus equations (5.44) and (5.45) give the vapour-pressure 

equations. 

Furthermore, at neghgible pressures, the free energy of the 
mixture is given by 


A 

kT 



= .YilogAi—.V,logV 


With equations (5.42) and (5.43) this gives 




where 


kT 


log giXi- Y,)-log Si-Sr log 4>r 

Yi log 6,-Sr log ir Prolog Yx'- -log Y/. 

4: log(Yi-gY,):-(4r- l)log(Yi-rY,)! 


and 


= -k(Yx-gY,)j(l-K)log 




kT 


(l-^)(l-/c) 


+ 


n(l_K)log 




(l-rj8)(l-rK) 


(5.46) 


(5.47) 


When the energy of mixing is zero, ^ is unity and A, vanishes. 
It mav therefore be regarded as the part of the free enerp- ue 
to the energy- of mixing. For the free energies of mixmg we 

obtain contributions 
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and 

kT 




+{l+/f)log 


2(i3+K) 

(l+)3)(l-fK) 


(5.49) 


These are the values already deduced by Orr in his treatment of 
this particular case; that equation (5.49) is identical with equa¬ 
tion (5.26) can be seen by substituting the value of X* in equation 
(5.24). 

The value of X on the assumption of the condition of quasi- 
chemical equilibrium is given in equation (5.20). Thus the heat 
of mixing is given by 

_ 1 _ „2 

A A W W A rV A V ^ ^ » 


AE = AH = Xw 


2(1+)3) 




(5.50) 


The entropy of mixing is given by 


TAS = AE~AA, 


(5.51) 


where AE is given by equation (5.50) and AA is the sum of AA 
given by equation (5.48) and of AA^ given by equation (5.49). 


Effect of non-zero energy of mixing on agreement 

with experiment 

5.11. The values of the entropy of mixing calculated in this 
way can be used to examine whether taking into account the 
non-zero energy of mixing improves the agreement between the 
experimental and theoretical determinations.f Prom the ex¬ 
pression for AG the partial potentials can be determined for 

vanous values of z and of t«. Typical results are shown in Fig. 10. 

From this diagram it follows that values of w greater than about 
500 calories correspond to highly insoluble systems; further¬ 
more, large negative values of w are unlikely for hydrocarbon 
systems. The best agreement with the experimental results, for 
reasonable values of z, is shown with w about 200 calories. The 
curves in Fig. 10 imply that in comparing the experimental and 
theoretical determinations of the entropy of mixing a value of 
w not greater than about 500 calories should be assumed. In 

t Orr, W. J. C., Trans, Faraday Soc, 40 (1944), 320. 

H 
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fact, such a value Tvill show the maximum difference which a 
Qon-zero energy of mixing can introduce. 

Making such calculations from equations (5.48), (5.49), and 
(5.50) TAiSi can he determined as a function of the volume frac¬ 
tion of the rubber. For volume fractions of rubber less than 



Volume fraction of rubber 


-1 


Fig 10. The free energj- of dilution of a polrmer-solvent 
«^-stem. The broken curve gives the smoothed exp^en^ 
r;sults for the natural rubber-benxene system. C^e (a) 
is for r = 5. If = 156 cal., curve (6) is for z - f 
cur%-e (c) is for z = 4. If = 358 cal., and curve (d) is for z - 4, 

ic = 536 cal- 

about 0-45 the value of T^S, calculated in this way is greate 
than that which is obtained when it is assumed the ene^ 
of mixing is zero. On the other hand, for r, greater than 0-4o tt 

value of TAS, is less than that calculated from 
theorv in which the energy* of mixing is 
is reckUed that the latter gave values of hig er ^ 
determined by experiment for dilute solutions and that „ 
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olutions rich in rubber, it will be seen that, over the greater part 
if the concentration range, the values of calculated from 
quations (5.48), (5.49), and (5.50) are not in quite so good agree- 
lent with the experimental values as those derived from the 
lementary theory which neglects the energy of mixing. The 



0*1 0-2 0*3 0-4 

Volume fraction of rubber 


Fig. 11 . Schematic diagram to show the effect of the 
energy of mixing on the entropy of dilution. Curve (a) 
is for a random mixture of molecules, curve (6) is for 
a mixtxire with a non-zero energy of mixing, and the 
broken curve shows the experimental values. 

;ate of affairs is shown schematically in Fig. 11, in which curves 
;) and (6) and the dotted curve respectively represent the results 
Dtained by the theories disregarding and taking accoimt of the 
3n-zero energy of mixing and the experimental results. It is 
ear that only between the points of intersection A and B of 
lese curves can the curve calculated from the theory which 
eludes the effect of a non-zero energy of mixing give better 
^eement with the experimental results. In particular, in the 
•gion of dilute solution where the discrepancy between the 
ementary theory and the experimental results is greatest, 
lowance for the energy change on mixing does notliing to 



5.12. It should also be noticed that the theoretical curve 
'hich is obtained when the non-zero energy of mixing is taken 
ito account is of the same shape as that for a mixture with zero 
nergy of mixin g. In particular, the theory which takes accoimt 
f a non-zero energy of mixing gives no indication of any point 


t/) 

< 


Volume fraction of rubber 

Fig. 12. The entropy of dilution of rubber-benzene 
mixtures in the region of dilute solution. Curve (a) 
shows the experimental results, curve (6) is the 
theoretical cur\'e for a random mixture, and curve (c) 
is for the energy of Tnixing (u?) equal to 536 cal.; the 
theoretical curves are calculated for s = 6. 

)f inflexion in the curves which are obtained by plotting 
is a function of the volume fraction of the polymer (tv) in the 
leighbourhood of r, = 0-25. It has already been seen that the 
point of inflexion at this point in the experimental curves arises 
[i-om the variation in the heat or energy of dilution deduced from 
sxperiment (Fig. 9). That, theoretically, a non-zero heat of 
movin g does not have such an effect would appear to cast doubt 
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on the validity of the variation of the heat of dilution as used by 
Gee (and shown in Fig. 9) and suggests that the shape of the curve 
for the entropy of dilution, deduced from these measurements, 
must be viewed with some reserve. In view of the wide scatter 
of the experimental points for the heat of mixing it is unlikely 
that a firm conclusion could be reached without more precise 
measurements of this quantity. In the meantime, however, it 
provides another reason, additional to those given in §§ 4.9 and 
4.10, why greater weight should be given to the agreement 
between the measured and calculated free energies for polymer- 
solvent systems than to the apparent divergence, for volume 
fractions of polymer less than about 0*25, between the measured 
and calculated values of the entropy. 

The introduction of a non-zero energy of mixing gives no 
better agreement with the experimental data than is given by 
the simple theory of a random mixture which takes account 
simply of the differences in size and shape of the molecules of 
the constituent species. Thus the heat or energy of mixing in 
itself cannot provide an explanation of the shape of the entropy 
curve for dilute solutions. It should, however, be noted that the 
analysis which has been given to include the effect of a non-zero 
heat or energy of mixing assumes that the configurational pro¬ 
perties of the polymer molecules are essentially the same for all 
volume fractions of polymer. This point will be discussed in 
§§ 6.3-6.5 below. 


5.13, Another important conclusion can be drawn from these 
calculations. Even with the maximum value of the energy of 
mixin g which is compatible with the existence of a solution, the 
maximum deviation of calculated from equations (5.48), 
(5.49), and (5.50), from the value which is predicted by the 
theory which assumes completely random mixing is not more 
than about 25 per cent. It follows that for a system in which the 
polymer is just soluble in the solvent the differences in the 


calculated values of the entropy are not likely to be more than 
about 20 per cent.; thus, the assumption of random mixing does 


not introduce a very serious error, particularly in comparison 
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tvith other possible errors which are inherent in the nature of the 
physical model which is used in the theory- of liquids. This 
inalysis provides a justification of the original assumption that 
n a first approximation it was reasonable to consider a system 
frith zero energy of mixing. It also indicates that the greater 
part of the divergence between the experimental results and the 
elementary theory which assumes random mixing must be due 
to some factor other than the non-zero energy or heat of mixing, 
lo a consideration of some of these factors and of the errors 


svhich are imphcit in the physical model of a hquid we must now 
:um in the closing chapter of this monograph. 
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CRITICISMS AND EXTENSIONS OF THE THEORY 

6.1- It has been seen that the introduction of a non-zei 
energy of mixing in the elementary theory of polymer-liqui 
mixtures developed in Chapter HI is unable to account for tl 
differences between this theory and the experimental result 
We shall now consider some other possible causes of the di 
crepancy between the theory of solutions of high polymers ar 
the experimental results. This examination will also serve 
indicate those aspects of the behaviour of solutions of polyme 
on which further research is required. Implicit in the applicatic 
of the theory to a particular system such as the rubber-benzei 
mixture is an assumption that the benzene molecule and tJ 
isoprene unit are sufficiently nearly the same size and shape th 
each can be allocated to one site of the quasi-crystalline arr? 
without serious error. Another source of error may be due to t] 
tendency, at least in very dilute solution, towards cluster form 
tion due to setf-coiLing of the rubber molecules. 

In this chapter we will first examine the effects of these t\ 
factors. Following this an indication of the effect of the appro3 
mate methods which have been used to construct the partiti* 
function will be considered. Finally, certain implications of t] 
physical model on which the theory is based will be considere 

The effect of allocating one site each to a benzene 

molecule and to an isoprene unit 

6.2. It has been noted that implicit in the application 
formulae such as equation (3.47) to the rubber-benzene syste 
is the assumption that the benzene molecule and the isopre: 
unit are approximately the same size and shape, so that ea 
can be allocated to one site without introducing any serio 
error. Bemalf has shown that for two sorts of spherical mol 
cules to pack together it is necessary that their molecular volum 
do not differ by a factor of more than two to one. Althouj 

t Quoted by Fowler and Guggenheim, Statistical Thermodynamics^ Ca 
bridge, 1939, § 814 fji. 
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Bernal's calculation is for spherical molecules, it probably mdi- 
cates approximately the volume ratio of molecules which depart 
not too greatly from a sphere and yet still pack together. In 
particular if the molecules of one kind are oblate spheroids whose 
shorter axis is the same as the diameter of the molecules of the 
other kind which are spherical, they -niU probably pack together 
if the other axes of the oblate spheroid do not differ from the 
shorter axis by a factor of more than about 1-4 to 1. In the 
course of another investigation Orry examined the available 
data, from X-ray analysis and a knowledge of the extent of the 
van der Waals field, on the sizes and shapes of the benzene mole¬ 
cule and the isoprene unit. His conclusions were that the benzene 
molecule in the hquid phase was roughly an oblate spheroid in 
shape having major axes of about 6-0 to 6-3 A and a minor axis 
of about 4-0 to 4-3 A with an average coordination of between 
6 and 8. and that the structure of rubber in the unextended 
state is probably fairly adequately represented by a roughly 
tetrahedral coordination of linked sub-groups with the distance 
between the centres of mass of neighbouring sub-groups about 
4-4 A. The ratio of the major axis of the benzene molecule to 
this spacing lies between T36 and 1-43. These figures suggest 
that the allocation in our phj'sical model of one site each to an 
isoprene miit and a benzene molecule may possibly lead to errors 
which would cause a divergence between the experimental and 
the theoretical results. 

To get some idea of the effect which this difference has on the 
thermod\Tiamical parameters of the system, the following model 
was examined, j The system was supposed to consist of poly¬ 
mer molecules as before and of disk-shaped molecules which 
are such that when one of them occupies a site, it precludes the 
occupation, by any other molecule, of those first shell sites which 
he in one of the planes (through the occupied site) of maximum 
abundance of closest neighbour sites. For instance, for simple 
lattices this means that each disk-shaped molecule occupies a 
site and four of its closest neighbours, while the other z—4 

t Orr, W. J. C., Trails. Faraday Soc. 43 (1947), 12. 

i Miller, A. R., unpublished results. 
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closest neighbour sites can be occupied b 3 ' other molecules. For 
various small values of r Dethe’s method was used to derive the 
partial differential equations which express the equilibrium 
of the system; these were then integrated to determine the 
corresponding combinatory factor; from these values a general 
expression could be written down (essentially a mathematical 
extrapolation) for any value of r. 

The formulae obtained conform to the general formula (§ 3.25) 
for mixtures of any number of polymer molecules of different 
sizes (equation (3.68) above), with each disk-shaped molecule 
regarded as a molecule of five segments corresponding to the 
number of sites which it has been assumed to occup 3 ^ At first 
sight this might appear to be somewhat surprising, but further 
consideration shows that the model of a disk-shaped molecule 
which has been used is really indistinguishable from a four- 
branched molecule which occupies a site and four of its closest 
neighbours with one branch going to each of these closest 
neighbours. This is only a special case of a branched chain 
molecule for which equation (3.68) is known to be valid, and 
consequently on the ph^’^sical model adopted for the disk¬ 
shaped molecules the general formula can be used with the 
appropriate values of the parameters. 

When the entropy of mixing is calculated for this case it is 
found to differ little from that calculated on the assumption that 
a polymer segment and a solvent molecule each occupy one site. 
In particular, it is no more successful than the elementary theory 
in reproducing the shape of the experimental curve in the region 
of dilute solutions. This can be made evident without difficulty. 
For the free energy of mixing of the polymer and the disk¬ 
shaped molecules the equation 


R 



( 6 . 1 ) 


is obtained, while the corresponding formula for a mono 

polymer system can be derived by rewriting equation (4.5) i 
terms of giving 

A(?i 


er- 
in 


R 



r->') ■ 


(6.2) 
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For the cases of interest r is very large, in particular r > d. 
Inspection of equations (6.1) and (6.2) shows that in this case 
and for dilute solutions calculated from equation (6.1) 
differs only insignificantly from calculated from equation 
(6.2). It must be concluded that the difference in size and shape 
between the disk-shaped molecules and the segments of the poly¬ 
mer affects the calculated thermodynamical functions hardly 
at all: a fortiori, the actual differences between a benzene mole¬ 
cule and an isoprene unit will be still less responsible for the 
differences between the elementary theory and experiment which 
have been noted in Chapter IV. 


Self-coiling of polymer molecules 

6.3, In polymer-liquid mixtures of moderate and high con¬ 
centration in polymer the individual polymer molecules inter¬ 
twine with one another so extensivelv that, even on a molecular 
scale, the availability of polymer segments is effectively unif orm. 
On the other hand, at very high dilutions the polymer chains 
separate out from each other and tend to form clusters by coUing 
up. Two lines of evidence point to such cluster formation. By 
studying dielectric relaxation in mixturest it has been shown 
that whereas normal mixtures show only a single loss maximum, 
mixtures which are near the point of separation into two phases 
show two maxima which mu:-t be due to considerable hetero¬ 
geneity on a molecular scale. It has also been shown! that for 
rubber-liquid mix tures falls rapidly with rv, that is, as 

more liquid is imbibed. This can only mean that few hquid- 
rubber contacts are made by molecules of liquid which enter the 
mixture near the saturation limit. 

The analvsis of mixtures containing dimer and trimer mole- 
cules using Bethe’s method allo^vs for the flexibility of these 
molecules. The formal mathematical extrapolation involved in 
the deduction of a formula appropriate to mixtures containing 
poh-mer molecules could be expected to yield correct results 
provided that no new effects are introduced when molecules of 
high polymers are considered. Furthermore, it was also indicated 

t SchaUamach, .4., Trang. Faraday Soc. 42 X (1946), 180. 

+ Qee, G. and Orr, W. J. C., ibid. 42 (1946), 507. 
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that provided the steric properties of the polymer molecules are 
the same in solution as in the pure species the additional con¬ 
figurations due to the flexibihty of the polymer molecules 
contribute nothing to the configurational free energy of mixing. 
With very long chain molecules, however, and in verj' dilute 
solutions, the flexibility of the poh’mer molecules may well lead 
to effects which are inhibited in more concentrated solutions. 


Such an effect is the seif-coiling of the polj'mer molecules, and 
it is not taken into account in the deduction of equation (3.47). 
Likewise in the derivation of the generalized combinatory factor 
for a completely random mixture by either statistical or kinetic 
methods, effects due to the self-coiling of the polymer molecules 
are disregarded in the assumption of complete randomness in the 
neighbourhood of a site occupied in a specified way. With very 
long polymer chains a segment of a polymer molecule may, as a 
result of the bending back of the molecule on itself, be competing 
for a site with a segment of the same molecule which has already 
been placed on the array of sites. The interpretation of complete 
randomness which has been adopted takes no account of such 


blocked configurations. This effect can, of course, arise only for 
long chains in which self-coiling occurs. It therefore appears 
that in very dilute solutions of high polymers aggregation into 
clusters occurs as a result of the self-coiling of the polymer mole¬ 
cules. The statistical treatments which have been given make 
no allowance for such effects. It therefore appears that the 
divergence between the theoretical and the experimental results 
in very dilute solutions is probably due to the coiling of the 
poljTner molecules with the consequent introduction of self- 
interfering or blocking configurations. 


6.4, On the basis of such a clustering of polymer molecules due 
to their coiling two attempts have been made to give a theory of 
very dilute polymer solutions. In one of these,! the region of 
space which is defined by a coiled polymer molecule is defined 
as the volume pervaded by a polymer molecule. A ‘swelling’ 
factor, in terms of which the thermodynamical quantities are 

t Flory, P. J., J. Chem. Phys. 13 (1945), 453. 
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measured, is defined as the ratio of the volume pervaded by the 
polymer molecule to its actual molecular volume; for flexible 
long chain molecules in solution this factor is very large. It is 
further assumed that, in very dilute solutions, those sites 'within 
the volume pervaded by the polymer molecule which are not 
occupied by segments of it are vacant. At moderate concentra¬ 
tions there is complete mtertwiniug of the pohuner molecules 
and every site is occupied either by a polymer segment or by a 
solvent molecule. It would appear that on such a model there 
must, at least for some concentrations, be a volume change on 
mixing which is contrary to the measurements which have been 
madej on the rubber-benzene system over the entire range of 
concentration. Consequently, this treatment is open to serious 
objection, and it appears that a more detailed analysis is neces- 
sars'. 

6.5. Orrj also considered the coiling of polymer molecules. 
As a first step in the analysis of the problem he examined the 
hmit of infinite dilution. In this limit the problem is that of 
determining the average configuration of a set of effectively non¬ 
overlapping flexible polymer chains in thermal equilibrium in a 
solvent phase. Ultimately the problem is to determine the 
statistical average configuration of a single polymer molecule in 
a given force field environment. Orr derived a formula for the 
heat of mixing which cont ains terms corresponding to the heat 
of solution of a hypothetical non-interacting chain and to the 
energ\- change in the transition from this hypothetical state to 
the equilibrium coiled-up configuration. Taking account of 
closest neighbour interactions and steric self-blocking effects, 
he computed the average configuration of polymer molecules; 
this calculation was carried out accurately by direct coimtmg 
for short chains but can only be done approximately, and in 
numerical terms, for long chains. Even so, this appears to be a 
more pro mis ing approach to the problem, for the following 

reason. 

t Posnjak. G., KoU.-Beihefte 3 (1912), 417; Gee, G., and Treloar, L. R. G., 
Tran^. Faraday Soc. 38 (1942), 158. 

i Orr, W. J. C., ibid. 43 (1947), 12. 
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Formally, the configurational term for a very dilute solution, 
in which the polymer molecules separate out individually and 
coil up, can be split up into a term to allow for the arrangement 
of the individual coiled polymer molecules and another term to 
give the number of configurations of each polymer on its cluster 
of sites on the quasi-crystalline array. It is the consideration of 
this latter problem that is the essence of this attempt to consider 
the limit of infini te dilution. Further progress appears to depend 
on its elucidation. 

6.6. It has been noted in § 5.11 that taking accoimt of a non¬ 
zero energy of mixing is unable to account for the point of in¬ 
flexion which is obtained in the experimental curves in which 
ASJvl is plotted as a function of the volume fraction of polymer. 
The theoretical cxirves which are obtained in this case are, in 
fact, of the same shape as those obtained for a random mixture. 
Thus ostensibly there is no reason on theoretical grounds why 
the heat or energy of mixing shoxild vary in such a way that it 
introduces a point of inflexion into the curve of plotted 

as a function of in the neighbourhood of = 0-25. It will be 
recalled, however, that the mathematical analysis given in 
Chapter V for a mixture with non-zero energy of mixing assumes 
that the configurational properties of the polymer molecules are 
independent of their concentration in the solution. In particular, 
no account is taken of the fact that it is likelv that in dilute 
solutions the polymer molecules cluster to form aggregates. 
The compaxison of the experimental and the theoretical results 
therefore suggests that the self-coiling of the polymer molecules 
in dilute solution, leading to their aggregation, affects the 
variation of the heat of dilution with the volume fraction of 
the polymer. This would account for the agreement between 
the free energy curves and the divergence, for dilute solutions, 
between the entropy curves. Thus it may be concluded that it 
is through the effect of the self-coiling of the polymer molecules 
on the variation of the heat of dilution that the divergence 
between the theoretical and experimental entropy curves arises 
for dilute solutions. Further research is therefore required to 
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examine the effects of the self-coiling of the polvnier molecules 
and, in particular, to examine how this affects the variation of 
the heat of dilution. 


The quasi-crystalline model 

6.7, It must be emphasized that the physical model of a 
liquid which has been used, and which is called a quasi-crvstalline 
model, is one in which a high degree of local order, without 
necessarily any long-range order, is envisaged. Furthermore, 
the statistical methods which have been used in Chapters III 
and V consider the detailed occupation either of a pair or of a 
small group of neighbouring sites. In no case does the argument 
depend upon the existence of long-range order. On the basis of 
the absence of long-range order in a liquid there appears to be 
little justification for any suggestion that the statistical methods 
which have been used, although valid for solid solutions, are of 
more doubtful applicability to liquid mixtures. 

These statistical methods do, however, assume that there is a 
constant coordination number which describes the local order of 


the quasi-crystalline array. How good an approximation to the 
partition function of the liquid mixture is provided by these 
statistical methods will depend on the extent of the fluctuations 
about the average coordination number. It is an essential 
assumption of the physical model of a hquid which has been 
adopted here, and which is currently accepted, that the fluctua¬ 
tions about the average coordination number, except perhaps in 
a negligible fraction of centres in the liquid, are far from able to 
disturb the regularity of the geometrical relationship between 
each molecule and its closest neighbours, that is, the short- 
range order. That this assumption is justified in a wide range of 
liquids has been verified by X-ray analysis.! 

6.8. Apart from this, the assumption of a constant average 
coordination number may introduce difficulties in the apphea- 


t Friedrich, W., Phys. Zeit. 14 (1913), 397; Debye, P., and Soherrer,_^P_., 
Gott Sach 16 (1916); Debve, P., and Meoke, H., Phys. Zai. 31 (1930), ; 

Ergeb. d. Tech, RofUgenk. H (1931); Bragg, W, L., The Crystalline State A 
General Survey (G. BeU and Sons, 1933), p- 195; Gingrich, N. S., Bev. Mod. 
Phys. 15 (1943), 90; Finbak, C., Avh. norske Vidensk. Akad. 10 (194o). 
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tion of the theory to particular systems, due to the fact that the 
pure liquids may have different coordination numbers. For 
rubber and benzene it has been shownf that the coordination 
numbers are probably a little greater than 4 and between 6 and 
8 respectively. The assumption of a constant coordination num¬ 
ber for the mixture may therefore be unduly restrictive in par¬ 
ticular applications. However, in view of the small dependence 
on the coordination number of the calculated thermodynamical 
functions, as shown by the curves in Fig. 7, this would appear 
not to exercise any great influence. 


Statistical methods 

6.9. It will be recalled that it has been shown that Bethe’s 


method agrees with the power series expansion in terms of 
moments of the configurational energy as far as the limit to 
which these have been calculated (§ 3.1). Thus while the Bragg- 
Williams approximation corresponds to the first moment of 
the configurational energy, the equations which are obtained 


using either Bethe’s method or the condition of quasi-chemical 


equilibrium correspond to an expansion as far as the third 


moment. 


6.10. These methods are, though, approximate. It has been 
shownj that the computation of the partition function for one- 
and two-dimensional arrays can be reduced to an eigenvalue 
problem. The method, which involves most elegant mathe¬ 
matics, appears to be particularly suitable for exa minin g those 
co-operative phenomena which exhibit transition points, such 
as order-disorder phenomena in alloys. It has the advantage 
that it can be computed exactly, but so far it has not proved 
possible to extend such analyses to a three-dimensional model. 
Furthermore, it assumes complete long-range order as in a 
crystal and consequently the method appears to be inapplicable 

t Orr. W. J. C., Trans, Faraday Soc, 43 (1947), 12, Appendix m; see § 6.2 
above. 

t Kramers, H. A., and Wannier, G. H., Phys. Rev. 60 (1941), 252, 263; 

Onsager, L., ibid. 65 (1944), 117; Wannier, G. H., Rev, Mod. Phys. 17 (1945), 

50; Onsager, L., Report on Conference of Phyeical Society at Cambridge, 1946. 
2 (1947), 137. 
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to liquid mixtures. This latter restriction can, however, be re¬ 
moved,! although in the apphcation of the method by which 

this is done to a three-dimensional model the results can onlv 

% 

be computed numerically. 

In this connexion two conclusions emerge from the considera¬ 
tions which have been examined in other sections of this mono¬ 
graph . The first is that mathematical refinement of the statistical 
calculation of the partition function is not the immediate need. J 
The second is that closer agreement with experiment is to be 
sought rather in constructing a physical model which provides a 
better representation of a liquid. 


Conclusion 


6.11. The behaviour of solutions of high polymers is markedly 
different from that of solutions of low molecular substances. 
Solutions of high polymers show large divergences from the laws 
of ideal solutions. The theorv which has been developed is based 


on a quasi-crystalline model of a Hquid, in which there is a high 
decrree of local order in the geometrical relations of a molecule 
and its near neighbours, but a complete absence of that long- 
order which is characteristic of a crystal. The local order 
is defined by an average coordination number, and fluctuations 
about it are neglected. This picture accords with the available 
X-rav evidence on the structure of hquids. The existence of such 
local'order and no long-range order is aU that is required for the 
apphcation of the statistical methods which have been used in 
developing the theory. Such a physical model cannot, though, 
claim to be precise and should probably be regarded as an 

approximation which is adequate for many purposes. 

The statistical theory of random mixtures which has been 
developed in this monograph agrees, at least semi-quantitatively, 
with the experimental data. It accounts for the greater part of 
the differences between the behaviour of solutions of high poly¬ 
mers and of ideal solutions. It establishes that the greater part 
of this difference in behaviour is due to the difference in size and 


hape of the molecules of the constituent species, 

t Domb, C., unpublished results. 

i Cf. Kirkwood. J. G.. J. Chem. Phys. 8 (1940), 6i3. 
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The agreement between the experimental and theoretical 
free-energy curves is very good. However, when the entropy 
is calculated, using measured values of the heat of dilution, it is 
found that the variation of the latter quantity is such that a 
point of inflexion is introduced, for high dilution, into the 
experimental curves in which ASJvj is plotted as a function of 
the volume fraction of the polymer. This is the point of least 
satisfactory agreement between the experimental results and 
the theory. In particular, introduction of a non-zero energy 
(or heat) of mixing into the theory fails to reproduce tliis effect. 
It has also been shown that this divergence is not due to the 
assumption that a solute molecule and a polymer segment are 
interchangeable and each occupies one site of the quasi-crj stal- 
line array, nor to the assumption of an average coordination 
number for the mixture. It should be pointed out, in this con¬ 
nexion, that the scatter of the experimental values of the heat of 
dilution is rather large, so that more precise measurements of 
this quantity would be desirable, before firm conclusions can 
be reached on this point. In the meantime greater weight 
should therefore be attached to the agreement between the 
experimental and theoretical free energy curv^es than to the 
divergence, in dilute solution, of the entropy curves. It is 
pointed out that there is evidence that in very dilute solutions 
the polymer molecules tend to coil up to form clusters. No such 
effect is taken into account in developing the theory which 
assumes that the configurational properties of the polymer 
molecules are independent of their concentration in the mixture. 
Thus it appears reasonable to conclude that the disagreement 
between the experimental results and the theory in the region of 
dilute solution is probably due to clustering of the polymer 
molecules, and, in particular, to the effect W'hich such fl-ggrega- 
tion of the polymer molecules has on the variation of the heat 

of dilution. Examination of this point is a subject for further 
research. 


That such good agreement exists between theory and experi¬ 
ment as has been shown in the preceding pages is evidence of 
the ‘striking advances in the mathematical treatment of the 

6129 

1 
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behaTiour of long-chain molecnlesT which have been made 
during the last decade. Further advances and more precise 
agreement are probably to be sought in a physical model which 
represents an actual liquid more faithfully than does the quasi- 

crystalline model, rather than in refinements of the mathematical 

% 

and statistical techniques. 

t Xature, 159 (1947), 734. 
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co-operative, 7, 28, 111. 

liquid as, 7. 
state of, 3. 

Bethe’s method, see Local configurations. 

Canonical ensemble, 4. 
grand, 6. 

Coiling of polymer molecules, 106, 113. 

effect of, on heat of dilution, 109. 

Combinatory factor: 

approximate methods of evaluation, see under particular methods. 
chemical kinetics, use of, 85, 88. 
defined, 6, 13, 16. 

difficulty of evaluation, 22 et seq., 27. 

efiect of configurational properties of molecules on, 43, 47, 58. 

evaluation equivalent to determination of partition function, 6, 22, 28. 

for mixtures containing dimer molecules, 39. 

for mixtures containing polymer molecules, 46, 57, 65, 92, 94. 

for mixtures containing trimer molecules, 43. 

kinetic derivation of, 59 et seq. 

partial differential equations for, 37, 56, 82. 

quasi-chemical equilibrium, condition of, 83, 88, 91. 

statistical derivation of, 55 et seq. 

Condensed system, 8. 

Configurations, direct enumeration of, 48 et seq. 
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configurational, 4, 5, 9, 10, 21. 
moments of, 28. 
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effect on agreement with experiment, 97, 113. 
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of concentration, 74. 
of dilution, 74, 78, 99. 
of mixing, 73 et seq., 97. 
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Free energy (coni,); 

of dilution, 63, 66, 75, 98. 
of mixin g, 14, 44, 63. 
for polj-mer solution, 46, 66, 88, 93, 96. 

independent of flexibility of molecules, 47, 51, 
partial molar, 63, 66, 74, 75, 98. 
of perfect solution, 13, 
power series expansion for, 28. 

Freezing-point, depression of, 64. 

Heat of dilution, 73, 80, 109, 

effect of variation of entropy, 80, 101. 
measurement of, 77, 80, 101. 

Henry’s law, 3. 

Interactions, closest neighbour; 

of a polymer molecule, 16, 19 et seq., 52. 
of various kinds, 21. 

Isoprene, 18, 103, 104. 

Kinetics, chemical, 83 et seq. 

Liquid model, 7, 18. 

isotropic harmonic oscillator, 11. 
quasi-ciystalline, 8, 18, 52, 110, 112. 

for polymer molecules, 18 et seq., 103, 110. 
smoothed potential, 10. 

Local configurations, method of, 28. 
detailed description of, 29 et seq. 

Mixtures, liquid: 

perfect or ideal, 1, 2, 13. 

Raoult’s law for, 2, 15. 
regular, 11, 15. 

Occupation of pairs of sites: 
probability, 54. 
relative firequencies, 53. 
used in method of local configurations, 30. 

Order: 

at a particular site, 29. 
local, 8, 110, 112. 
long-range, 9, 111, 112. 

Osmotic pressure, 64, 68. 
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Partition function (coni .): 
grand, 6, 

of nuxtures containing dimer molecules, 33. 
of mixtures containing poIjTner molecules, 17, 21, 65. 
of mixtures containing trimer molecules, 39. 
of perfect solution, 13. 
of regular solution, 11 et seq., 31. 
of assembly of molecules, 9. 
of pure liquid, 13. 

Polystyrene-solvent systems, experimental, 69, see Vapour pressure. 

Raoult’s law, 2, 15, 66, 68, 70, 73, 81. 

Rubber-benzene systems, 72 et seq., 97 et seq., 104. 

Swelling pressure, 64, 

Vapour pressure: 
equations, 15, 65 et seq., 96. 

comparison with experiment, 70 et seq. 
experimental determination, 64. 
of polystyrene-solvent systems, 70, 71. 
of rubber-benzene system, 72. 
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